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NATIONAL ADVISORY COMMITTEE FOR AESOIAUTICS 


TECHNICAL NOTE NO. 1115 


INVESTIGATION OF THE STABILITY OF THE LAMINAE 


BOUNDARY LAYER IN A COMPRESSIBLE FLUID 


By Lestor Lees and Chia Chiao Lin 


SUMMARY 


-i.n the present report the stability of two-dimensional laminar 
flows of a gas is investigated by the method of small perturbations. 
The chief emphasis is placed on the case of the laminar boundary layer. 


Part I of the present report deals with the general mathematical 
theory , The genera,! equations governing one normal mode of the 1 small 
velocity and temperature disturbances are derived and studied in great 
detail. It is found that for Reynolds numbers of the order of those 
encountered in most aerodynamic problems, the temperature disturbances 
have only a negligible effect on those particular velocity solutions 
which depend primarily on the viscosity coefficient ("viscous solu- 
tions"). Indeed, the latter are actually of the same form in the 
compressible fluid as in the incompressible fluid, at least to the first 
approximation. Because of this fact, the mathematical analysis is 
greatly simplified. The final equation determining the characteristic 
values of the stability problem depends on the "inviscid solutions" and 
the function of Tietjens in a manner very similar to the case of the in- 
compressible fluid. The second viscosity coefficient and the coefficient 
of heat conductivity do not enter the problem; only the ordinary coeffi- 
cient of viscosity near the solid surface is involved. 


Part II deals with the limiting case of infinite Reynolds numbers . 

The study of energy relations is very much emphasized. It is shown that 
the disturbance will gain energy from the main flow if the gradient of 
the product of mean density end mean, vorticity near the solid surface 
has a sign opposite to that near the outer edge of the boundary layer. 

A general stability criterion has been obtained in terms of the 
gradient of the product of density and vorticity, analogous to the 
Rayleigh-Tollmien criterion for the case of an incompressible fluid. If 
this gradient vanishes for some value of the velocity ratio of the main flow 
exceeding 1-i/M (where M is the free stream Mach number), 
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then neutral and self -excited "subsonic" disturbances exist in the in- 
viscid fluid. (The subsonic disturbances die out rapidly with distance 
from the solid surface.) The conditions for the existence of other 
types of disturbance have not yet been established to this extent of 
exactness. A formula has been worked out to give the amplitude ratio of 
incoming and reflected sound waves. 

It is found in the present investigation that when the solid bound- - 
ary is heated, the boundary layer flow is destabilized through the change 
in the distribution of the product of density and vorticity, but stabi- 
lized through the increase of kinematic viscosity near the solid bound- 
ary . When the solid boundary is cooled, the situation is Just the 
reverse. The actual extent to which these two effects counteract each 
other can only be settled by actual computation or some approximate 
estimates of the minimum critical Reynolds number. This question will 
be investigated in a subsequent report. 

Part III deals with the stability of laminar flows in a perfect gas 
with the effect of viscosity included. The method for the numerical 
computation of the stability limit is outlined; detailed numerical cal- 
culations will be carried out in a subsequent report. 


INTRODUCTION 


In a recent paper (reference 1), one of the present authors has 
clarified the theory of the stability of two-dimensional parallel flows 
in a homogeneous viscous incompressible fluid. The experimental inves- 
tigations of H. L, Drydon, G. B. Schubauer, H. K. Skramstad (reference 
2) and H. W. Liepmann (reference 3) agree with the calculations made by 
Tollmien (reference 4) , Schlichtlng (reference 5) and those given in the 
paper quoted (reference l). Because of the increasing importance of 
phenomena of gas flow at high speeds, it seems natural that the investi- 
gation should be extended to cover the case of a gas, taking into ac- 
count the effects of compressibility and heat transfer. 

The interest in this problem is further enhanced by the fact that 
disturbances of finite amplitude in high-speed flows are known, to have 
the tendency of building themselves up into shock waves. It is there- 
fore possible that instability of high-speed laminar flows will lead to 
shock waves instead of turbulence. Although an instability theory in- 
volving only small disturbances would not be able to settle this point, 
it at least paves the way to such an investigation. 


This possibility was first pointed out to the authors by Doctor 
H . W. Liepmann. 


2 
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The present report is concerned with the stability of two- 
dimensional laminar flows of a perfect gas, subject to small disturb- 
ances. The chief emphasis is placed on the case of the boundary layer. 

It is hoped that the results may throw light on the general features of 
the relation of compression waves with the boundary layer, that it my be 
known how the . stability of a boundary layer is affected by the free 
stream velocity, and by the thermal conditions at the solid boundary. 

As in the incompressible case, only small two-dimensional wavy dis- 
turbances will be considered. Unfortunately, there is no rigorous' proof 
in the present case that these disturbances are more unstable than the 
three-dimensional ones. But the results of the incompressible case to- 
gether with some physical considerations seem to justify such a treatment, 
whicu naturally simplifies the mathematical analysis to a large extent. 

One essential difference of the present problem from the install 1~ 

. y problem in an incompressible fluid is the presence of an appreciable 
interchange of .mechanical and heat energies. Another is the fact that 
the flow velocity is of the same order of magnitude as the velocity of 
sound. The present investigations, however, reveal that the chief phys- 
i cal mechanism, is not changed. That is, the stability of two-dimensional 
parallel flows depends primarily on the distribution of angular momentum 
01 an element of the fluid, and on the effect of viscous forces, but not 

° n * lea 'k con Activity . The viscous forces influence the stability 
o the flow both in building up the disturbance by increasing the 
Reynolds shear ^stress and in destroying the disturbance by dissipation. 
(Cl. sec. 14 of reference 1.) In the present case, however, the "angu- 
lar momentum of a given v olume of the fluid depends upon the product of 
.4 g _ n ff i ty tjcity . Thus, the gr ad ient of t his product plays the 
same role as the curvature of the velocity profile (gradient of vortic- 
ity) in the incompressible case. Moreover, since the magnitude of 
viscosity varies with temperature across the flow, there is an uncer- 
tainty in defining a Reynolds number which will properly describe the 
stability characteristics under various conditions. It is concluded 
from the present investigations that the viscosity coefficient in the 

important . This tends to justify 
the process of Alien and Hitzberg (reference ' ~6)~ in estimating the criti- 
cal Reynolds number for the boundary layer of a compressible' fluid, so 
far as their treatment of the viscosity coefficient is concerned. They 
have, however, neglected the effect of the distribution of angular 
momentum in the fluid. 

It i found in the present investigation that when the solid bound- 
ary ls heated, the boundary layer flow i s de stabilized “through the 

of distribution of angul ar mome ntum, bu Fs^l'lize^^ough t he 
1 i^j , _ ease of kinematic viscosity n ear the solid boundary. Mien the 'soli d 


3 
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boundary is cooled, the situation is just the reverse . The actual ex- 
tent to which these two effects counteract each other can only he settled 
hy actual computation or some approximate estimates of tho minimum crit- 
ical Reynolds number. This question will he investigated in a subse- 
quent report . 

The transfer of energy between the mean flow and the disturbance is 
also somewhat different from that in the incompressible case, because 
the velocity of sound is finite . When a disturbance is being amplified, 
energy passes from, tho main flow to tho disturbance not only to supply 
the increase of energy of disturbance inside the boundary layer, but 
also to supply the energy carried out of the boundary layer by the disturb- 
ance. For a damped disturbance, the opposite is true. For supersonic 
flows , there is also the possibility of energy transfer for neutral dis- 
turbances, The energy is carried into the boundary layer by an incoming - 
wave and. out c.i: it by an outgoing wave. These two waves are generally 
present simultaneously, ana the situation may be described as a reflec- 
tion with a change of amplitude , 


From the behavior of the disturbances outside the boun dar y layer, 
they aro classified as "subsonic,” "sonic," and "supersonic” disturb- * 
ances, according to wnother tue phase velocity of the wavy disturbance, 
in the direction of the free stream and relative to an observer moving 


with the velocity of the free stream, 


local velocity 


7 o:i: 


sound. 


-s below, equal to, or above ti 


The method of analy Loj-b & xxi tne px*6 souls 1.8 vvi v y sil.rxiX^x* 

to that used in the incompressible case. Indeed, an attempt is made to 
establish results analogous to those obtained in that case. Small dis- 
turbances are considered, which are analyzed linearly into normal modes, 
periodic in the direction of the free stream, and these are then tree! ' 
separately. The normal modes may be damped, neutral, oi‘ self-excitod 
oscillations in time. TiVw> « <-oh-v~ — ----- *■> -- - - 


For a given condition, the main flow 




unotak 
ince be- 
fci.le 


if any one of these 

becomes so large that it can no longer be regarded as linear, 
present theory ceases to apply. But it may then be expected that turbu- 
lenco or shock wave will be precipitated by the nonlinear effect. 

Reference should be made to the work of Scklichting (reference 7) 
and Ku.chema.nn (reference 8). Schlichting was interested in the stabi- 
lizing or destabilizing effect of gravity and temperature gradient. But 
he neglected the interchange of mechanical and heat energies. In aero- 
dynamical problems, the authors are not particularly interested in the 
el x ect of gravity. However, in the general, mathematical investigation, 
the inadequacy in Schliclitiny’s fiuulaiaental equation will be discussed 
(sec. 2).. ICiicheman made only an "inviscid" investigation of the stabil- 
ity of the boundary layer (pt. II of this report) , but he neglected the 


4 
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gradient of temperature and the curvature of the velocity profile. These 
are serious omissions. Their significance will turn out in the process 
of the present investigation. 

This investigation, conducted at the California Institute of 
Technology, was sponsored by and conducted, with the financial assistance 
of the National Advisory Committee for Aeronautics. 

The authors are indebted to Professor Theodor von Karman for super- 
vising the work and to Professors Clark B. Millikan and E. S. Tsien for 
their interest and discussions. 
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Charac- 

teristic 

Dimensional quantities Dimensionless quantitie s measure 


Components of stress tensor 


(9) 

* -* *t 

T = T +T 

XX XX XX 

» 

T as T + T 

XX XX XX 

*o 

(10) 

* •** #t 

T * T + T 

xy xy xy 

I 

T = T + T 

xy xy xy 

Po 

(11) 

* --K- #t 

T = T + T 

yy yy yy 

» 

T * T 4* T 

yy yy yy 

_ -)f 

Po 

Density of the gas 



(12) 

* _->■ *» 

p = p + p 

/ v /v ia(x - ct) 

Ply) + r(y)e 

Po 

Pressure of 

the gas 



(13) 

* _* *» 

P * P + p 
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Temperature of the gas 



(UO 
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0 
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of viscosity of the 
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tij(y) + mi(y)e ic, ( x " ct ) 
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Dimens i onal q uantities D imens i onless quantities 


Charac- 

teristic 

measure 


Thermal conductivity 


( 18 ) 


-X- _* * * 

k = k + k 


= ni(y) +-^~k(y)° la ^ x " ct ^ c p Pio 


Wave number of the disturbance 


(19) 


* * 

a = 2n/\ 


a ss 2n/\ 


-i 


Phase velocity of the disturbance 


( 20 ) 


X) 


Specific heat at constant volume 


(21) C v 1 C v 

Specific heat at constant pressure 

(22) C p 7 Cv 


Gas constant per gram 

■Sr 

(23) E 7-1 C v 


Acceleration duo to 'gravity (in the negative direction of the y-axis) 


(24) g 



2 


U Q / 1 
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Dim e nsion less quantities 


Froude number 


(25) 


* 



■Reynolds number 


(26) 


R ~ Po u o-/^xo 


Mach number 



M = 



Prandtl number 


(28) 



ior the case of the boundary layer, the boundary- layer thick- 
ness 5 vill in general be taken to be the characteristic length; .for 
some purposes, the displacement thickness 8j_ will be used. A bar over 
a quantity denotes average value, a dash denotes fluctuation; and the 
subscript ( ) Q denotes free stream value in the case of the boundary 

j-ayer. The subscripts ^ and p denote the real and imaginary 
parts of a quantity, respectively. 


I - GENERAL THEORY 

1. The General Equations of Disturbance 


The general equations of disturbance for a perfect gas which is 
flowing parallel or nearly parallel to a given direction will now be de- 
rived. . As has been explained, only two-dimensional motions with two- 
dimensional disturbances will be considered. 


8 
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With the system of notation explained above, and with positive y-axis 
pointing vertically upward, the general equations for two-dimensional 
laotion of a perfect gas may ho written as follows: 


(a) Equations of motion, 


Su 

St* 


+ u 


* Su 


Sx' 


+ v 


* d u _ Y d T xx + ^ T xy 

* * ' V - N * * 

cy p ox Sy 


( 1 ) 


Sv* 


. Sv* Sv* 

+ U* + V* 


St* OX* 


Sy* 


1/ zz V g* 

p * ' v Sx* Sy* ^ 


( 2 ) 


(h) Equation of continuity, 


M JL 

St* + Sx* 


( p*u* ) + ( p*v* ) = 0 

Sy 


(3) 


(c) Equation of energy, 



(d) Equation of state, 

p* = p*R*T* (5) 


In these equations, exx, e^y, €yy are the components of the rate -of - 
strain tensor, and t^, f'xy, t yy are the components of the stress 
tensor. They are defined as follows: 


9 
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ha* 


•xx 


ox * ' 


* = 1_/%L+ hrf\ 

e xy 2 W c&V ’ 


*• 

e yy 


dy* 


( 6 ) 


* 

T xx 


-p" 


* 3u 
2ui — 


* 


_ _§./ * 

8x* + 3 


Hi) 


/du* 

V8x* 


frr* \ 
+ ay V 



-x/ 3u* 



Svf\ 

+ dxV 



■X- ^ 

P + 2Ui 



( 7 ) 


J 


The coefficients of viscosity n* and \i% and the coefficient of 
heat conductivity k'£ are essentially functions of temperature. Hence, 
there results a system of five differential equations for the five 
variables p*, T*, p*, u*, v*. 

Consider a motion which is slightly disturbed from a steady state. 
Then it is convenient to separate any quantity Q*(x*, y*, t x ) into a 
steady- state part Q* (x*, y*), and a small disturbance Q* f ( x* , y*, t*) 

Q*{x*, y*, t*) = Q*(x*, 'y* ) + Q*‘(x*, y*, t*) (8) 


By substituting expressions of the type (8) for each of the variables 
into (1) to (7), remembering that the steady- s.tate parts satisfy those 
equations by themselves, and, finally, neglecting terms quadratic in 
the small disturbance, the following system of equations of disturbance 
is arrived at: 


du* f „ , 8u* 


St" 


4- \ U 


* * 


jhr \ 


+ v ->' 1 j + 

hx* by*S 


u* 


<3u 


* t 


du* 


-f. v 

ox* 


V 


,/8t* x . c> 


* In 




dX 


dy* 


*L) - 


P* /CiT. 




c>.t„ 


xx , w xy 


P ^ dx[* c)y x ' 


(9) 
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***'.[.*'<** *t&*\ l 


+ !u* 

St* \ 


dx 


+ v + tl + V 


Sy*/ 


x - * x ** 
ov »dv 


Sx 


'X* X 3$* 


i * t * t 

1 1 S St 




77 


Sy 

*t , «* 

p /St^ Sr yy 


-* \ -n * •% * 

p \ Sx Sy 


p X \ Sx* Sy* 


(10) 


*F 

sp + _s_ 

x # 

ot Sx 


f ’ * -* **\ A / -•<•» _* _-x- *«\ 

v^P u + P u y + y^p v + p v • j = 0 

Sy* 


( 11 ) 


/ vr * * f 

n -*/ ST -* St 

Co + u 

V \ * VC 

\ St dx 


*t 

st St ** ST 


+ V — + U + V 

•v n 

S y Sx Sy 


\ 


+ c v p*'/ il + u* iL + ;* il) „ *-(i**L- + k*’i^ 

dt S- Sy / Sx* \ ax* Sx* 


\ 


S /-*ST*' 


Sy \ Sy Sy 


_# 

** St \ /-* ** _* *» _* *» 

+ k — + ( v T xx c xx + 21 xy e xy + T yy e yy 


* » _* *t _# *» _* 

+ ( T XX e XX + 2T xy e xy + T yy e yy 


( 12 ) 


* t *1 -x - 1 

p p r 


-* _-x- 




P P T 


(13) 


where 
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(Relations of this type, of course, do not hold for the disturbances.) 
By considering only a local region, say around x* = x* 0 , and introduc- 
ing the boundary- layer approximation, the flow is regarded as essen- 
tially parallel, with every mean quantity q' evaluated at x* - x* 0 . 

Thus, for parallel or nearly parallel flows, the differential equations 
of disturbance do not contain x* and t* explicitly and an attempt 
may be made to find solutions of the type 


* * t * -X- 

Q t (x. 9 j > • 


t ) 


(y * )e ia (x -c' t ) 


(IT) 


Indeed, every quantity will be reduced to a dimensionless form in 
accordance with the scheme listed (List of Symbols); for example, 


f (/) = u Q w(y ) , 


-v- I 


u 


= U- 


f (j)o ia<X - Ct) 


i,3) 


For the present, the additional restriction 
velocity temperature f*, and so forth, vili 

£ir^l_.dipiens and, linear! zed epiuatioris for 

distur bances) in a par allel or n early paralj e j. mai 


of unif ona free stream 
be retained* Then the 

th.e_ amplitude b of small 
n flou are as foil or n : 


a p 


ii(¥ 


x f I 

c)f + v epf 


lew? 
+ 






? h2 - m-x .... r 1 

+ - - at 9-f + im‘ T 

3 S L 


+ — • mxw* ’+ m* w* + |i£ (f 1 + iaf’cp) !>• 

R l ' J 


(19) 
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a“'p -f i(w - c)cp =i . .. 


t 


/ 


Pa.a, 

+ J2q>' ' - if ' - a2q> 


A " 


yM" F r ^ 


2a ^ ^ ( t , , ■] 

+ J. ( p + if L 

3 B l J 


<X f ftp ^ 

+ jfi im iW * + 2piCp + j. (n 2 - Hi) (<p* + If) j. 


( 20 ) 


i(w - c)r + p(cp 1 + if) + p'cp = 0 


(21) 


ap i(w - c)0 + T ’cp {- = -a( y - l)pT(cp f + if) 

V J 

y r 


^ Hi(® - a^Q) + (m.iT ) + j-i x 0 7 


T t 


J t J * 1 


Ea r 


f y « 1 ) ( 1 2 f | ! 

+ y 1 M 2 7 m^w + 2[iiw (f + ia 2 cp) 4 


E 


l 


J 


( 22 ) 


it r 6 

P P + T 


(23) 


The following two dimensionless equations for mean quantities should 
also he noted: 


Equation of state 


p = pT 


(24) 
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Equation of static pressure gradient across the boundary layer. 


2 

i 7M 

p = p (25) 

F 2 


J.n^alx these equations, and in all subsequent equations, a dash denotes 
differentiation with respect to the dimensionless variable y, and 
should not be connected with the idea of a fluctuation. For example. 


f dpi d Hi dT 

= B ( 2 g) 

dy dT dy 


while the fluctuation mi is given by 


mi 



(27) 


Note that in (18) a characteristic velocity has been used as the 
reference variable. This stresses the role of the inertial forces. In- 
deed, the simu-Lbaneous comparison of the inertia forces with pressure 
and gravitation, as embodied in the Mach number 


_* 



* 


and the Froude number 
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makes it difficult to consider the limit of small inertial forces 
u 0 •—» o. For in such a limiting case, the phenomenon is essentially 
governed by pressure and gravitational forces, which both become in- 
finitely large compared with the reference inertial force. In aero- 
dynamical problems, however, this limiting case is not of importance. 

Mathematically speaking, a singularity is brought into (20) if the 
neighborhood of F = 0 is considered. Thus, it is possible to study 
only the case of a small Mach number. The limit of vanishing Mach num- 
ber and vanishing Froude number can be considered only when r = 0. 
These statements will become clearer after reading the detailed discus- 
sions in section 3. 


To study the case of an extremely large Mach number, on the other 
hand, ^it would be more convenient to use the stagnation or "rest" val- 
ues of pressure, density, and temperature as characteristic measures, 
rather than the free stream values. 


2. Analytical Nature of the System of Equations 
of Disturbance and Its Solutions 

The system of equations of disturbance (19) to (23) consists of 
five linear equations in the five variables f, 9 , n, r, 0, with p, 

P, T, w supposedly known from the steady-state solutions. Before " 
^-pply-^g this system to any definite problem, it is necessary to know 
clearly its analytical nature; for example, the number' of sets of lin*- 
early Independent solutions it possesses must be known. It is also 
desirable to know the general analytical nature of the solutions in 
the variable y and in the parameters M-, F2, E, a, and d. In all 
these discussions of analytical nature, both the variable y and the 
parameters will be regarded as complex. 

To settle these questions, it is convenient to choose a number of 
new variables 'Lx, . . ., 2>n and rewrite the system into the form 

. n 

dZi 

- / t (y) z j (i = 1, 2, . . ., n) (28) 

J=i 


where A ± j are known functions of y. Since (19), (20), and (22) 
involve the second derivatives of f, cp, and Q, it seems desirable 
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to choose the six dependent variables as f, <p, 0, f*, y* } 0*. j n 
this way, the equations can be set equivalent to six equations of the 
"type (28), if r and n are supposedly solved algebraically from (21) 
and (23) ; It is seen, therefore, that the system of five equations (19) 
to (23) is actually equivalent to six homogeneous linear differential 

equations of the first order, and there are six linearly independent 
solutions . 

a However, this choice of the dependent variables is not satisfac- 
tory. It leads at once to the suspicion that the solutions have 
singularities at the point where w = c. Fop, in solving for r from 
(cl) a singularity is introduced into the coefficients j ( y ) of the 

system (28). Phys.ically, the solutions cannot have such a singularity 
for real values of y. Hence, it is necessary that such a singularity 
De ; Qnl y apparent. Indeed, this can be shown to be true by a new choice 
of the dependent variables. 

To be more precise, let it be assumed that the known functions w, 

p, p, T are analytic functions of y and of the parameters M 2 , — . 

F 2 

These j. unctions may be regarded as independent of the Reynolds number 
\, when .the characteristic length l is properly chosen. This assump- 
tion is related to the boundary-layer approximation, and is therefore 
accurate up to the same order. For example, for the Blasius profile, it 

is accurate up to the order of (R5) 5 being the thickness of the 

boundary layer. < 

How choose the system of dependent variables 


Zi = f. 


Then, at once. 


Z2 = f ^ , Z3 = cp 


Zs = 0, z 6 = e» 


(29) 


dZiL 

dy 


= Z 5 


(30) 
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~ = z 6 (31) 

ay 


The equation (23) becomes 




(32) 


.L J. 


vhich makes it possible to eliminate the variable r without increa 
the order of the differential equations and without introducing? anv 
singularity. Then from ( 19)., (21), and (23), solve for z' 2 , Z' 3 ,"and 
’ e which, when reduced with the help of ‘ ( 30 ) to (32) are equations 

of the type ( 28 ) with A i ^(y) regul ar in both y and the parameters. 

The equation for Z 0 is 


Z 




i(* 


, i i-r 

: ) 2 ~~ Z/ 

Ip 4 



.33) 


but the other two equations are too lengthy to be written out explicitly. 
They are of the following general nature: 


dZ 2 aR 
!y Hi 



c)Zx + w'z.. 


+ tz 4 r + o(i) 


r T7 

dy 


ooE 
7 hi 


r 

-•lyp i(w 


i. 


c)Z 5 + 



- (7 



i(w 



!• + 0 ( 1 ) 
J 


(34) 


33) 


where 0(1) denotes a linear function of Z x> . . . , Z 6 which is of the 
order of unity in the parameter E and is regular in the parameter M ;/: . 
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Tlie differential equation for Z 4 must "be obtained from (20) in a 
slightly different manner. It is necessary first to eliminate <p" by 
using (33) and then solve for n ' = M 2 Z '. Tiiere is obtained 

4 




r 

'i 


1 + yM‘ ; 


io 2 , 

Qi 2 + 

B 3 




(w 


P 


o)T 


1 


x f - -- (-?- m 2 z 4 - £l z . V 
l Vp T / 


ia(w - c)pZ 3 +— 0(1)1 

B J 


(3 6 ) 


where 0(1) has the same general meaning as before. It is noted that the 
last step is the only division involved in this process of elimination. 
Thus, unless 


1 + - -((ig + 2|ii) 

E 3 



(37) 


(which is not possible for I R S » l), the system of differential equa- 
tions (30) to (36) is reg ular in y and in th e paramet ers. But since 
the regularity breaks 'down for infinite R and for R satisfying (37), 
any expansion of the solution as a power series in R must be in the 


form of Laurent series . In the parameters — - and M 2 , the coeffi- 

F‘~ 

dents are entire functions; in the parameters a and c, they are 
analytic in a region including the origin. 


From the general existence proof of the solutions of linear differ- 
ential equations by means of successive approximations, it is clear that 
these properties of the coefficients persist in the solutions. That is, 

there exists a fundamental system of six solutions Zj_(y; M 2 , l/F 2 , 

R, a, c) (i s 1, 2 , 3, 5, 6) which are analytic functions of 

y and of the parameters. 


Now consider a few limiting cases: (l) M 2 ~$>0, (2) l/F 2 — >0, 

(3) E > « . As discussed at the end of section 2, if M 2 — -$> 0 by 
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making the velocity &* — > 0, then F 2 and R approach zero at the 
same time. This is certainly not what the authors wish to discuss. Rath- 
er, they are thinking of the solutions ,Z i as expanded in power 

series of the parameters, say. 



and retaining only the zeroth order term as an approximation. This 

process is valid so long as M 2 is sufficiently small compared with 

unity , while hoth F 2 and 1? are of their usual magnitudes (namely, > 

much larger than unity). For convenience, the mathematical process 

M 2 — $> 0 will still be preserved. But this must not be confused with 

any physical requirement that u — > 0 or c* — > •<* . the latter 

o o 7 

being in contradiction with the equation of state. The limiting case 

R oo is an asymptotic approximation and will be dealt with more 
• carefully below. 

Case (1) M 2 — >0 . With the relation ( 24) in mind, the equations 
(19) to (23) become (ca » it/yM 2 ) 


f . . * ^ f # t 

ap<|i(w-c)f + w<p^* iam + -l f + a 2 (icp - 2f ) 

J R l 


, 2 H 2 - Hi 
+ - ■ .. . .. 

3 B 


•t + icp 




r 


R 


t 


mjv + miw + 


ul(f’ 


ia 2 cp)j> 


(39) 
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r 


i(, - c)(fl . -4 - If - lap] * ~ * if j 


C F ' E l ' j 3 • E 


+ “|"im 1 w' + 2 hi<p’ + *§■ (us - UiM 1 ?' + if) 1 

E 7 3 j 


(40) 


i(w - c)r + p(cp* + if) + p'cp = 0 


(41) 


r h 

ap < i(w - c)d + T'cp | = -<x(7"l)pT(cp' + if) 

J 


Be, 


• | ux(e" - a 2 0) + (h-i 1 ’ ' ) + Hi® 


(42) 


r o 

o =_ + — 
P T 


(43) 


The fact that equation (25) reduces to p' = 0 indicates that the grav- 
itational force is important in these problems only insofar as the 
buoyancy corresponding to density fluctuations is concerned and not in 
connection with the determination of mean density distribution. 

This set of equations (39) to (43) is different from that used by 
Sclilichting (reference 7); who neglected temperature variations but in- 
cluded density variations. In his case, (42) yields the condition of 
i nc ompre s s ibi 1 i ty 

qr + if = 0 (42a) 


and (4l) becomes 


i(w - c)r + p'cp = 0 
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He also made certain other minor reductions in (39) and (40). The jus- 
tificatior of Schlichting’s assumptions is not obvious. Also, his 
complete equation of disturbance (equation (11), p. 319, reference 7). 
has a singularity at the point where the phase velocity is equal to the 
mean velocity of the flow. This gives rise to multi pie -valued solutions, 
to which it is difficult to assign a proper physical interpretation. 

Reduction to the equa ti on o f Orr and_ Somerf eld for an homogeneous 
incompressible fluid .- This simple case is obtained from the limiting 
case of zero Mach number with the additional requirements that the mean 
pressure, mean temperature, and mean density are constants. These con- 
ditions can hold only in the case where there is no conduction of heat 
across the boundaries. Otherwise, there must be a finite - temperature 
gradient at the boundary. When p, p, T are constants, indeed 
p = p = T = 1; then the equations ( 4.1) to (43) give 


i(w - c)r + (9’+ if) = 0 (44) 


r + 6 = 0 ( 45 ) 


ia(v - c)8 - — p 10 (G - af’C) (46) 


Multiplying (46) by 0, adding the corresponding complex conjugate, and 
integrating between the boundaries along the real axis of the y- plane 
gives 


ya 


e) a 


0 , |' dy - 2(oA + cI!o 0 Cj /p-io) / C j ' dy 


yi 


yi 


iil tlio bounciary 
n o 1 ;u t i on o i’ ( 4 6 ) 


onditioiiG are a* ■ 
ie not identically 


at Doth houndarlo;.; • 
:*o only alien 


the 


c i < - aHio/d^ 


(47) 
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Thus, if the main interest is in the limit of stability (c^ changing 
sign), only the solution 6=0 need he considered. Then r = 0 by 
(45), and the equation (44) reduces to 


cp' + if = 0 


(44a) 


The equations (19) and (20) then become 


1 


ia(w - c)f + aw'cp = -iaco + — (f - a£f) 

E 


> (48) 


ia 2 (w - c)cp = -cd' + — (cp" - a 2 cp) 

E 


Eliminating cd from these equations and then substituting f from 
(44a), the equation of Orr and Sommerfeld is obtained. 

Case (2) ff 2 - -> cjo . In this case, all the equations (19) to (25) 
remain unaltered, except that the terns in l/F 2 should be dropped from 
( 20 ) and (25). This is the case which will be discussed more in detail. 

Case ( 3 ) E — > o> . In this case, the equations (19) to (23) become 


P f i(w 


c)f + w'cp'’ - - i?- 

J y M 2 


(49) 


ia 2 (w - cOcpl = 


yW 


r_ 

F s 


(50) 


i(w - c)r + p(cp' + if) + p'cp = 0 


(51) 
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c) 0 + T'cp 




( 52 ) 



( 53 ) 


The equations (24), (25) remain unchanged. It is to be noticed that the 
orders of the differential equations are reduced. This is consistent 
with the fact that the solutions have an essential singularity at 

E = co. 

After the elimination of f , it, r, 9, the final differential 
equation for cp reads 


d f(w - c)cp' - w'cjJ] of ; ( w - _ 1 

dy T - M s (v - cWJ T ’ F 


( 54 ) 


where L (cp). is a linear expression in cp involving cp* and cp. 

The boundary-va lu e problems, - For a given physical problem, there 
are usually associated certain boundary conditions on the disturbance. 
For example, for flow between fixed parallel plates, the velocity dis- 
turbances must vanish at these plates. Also, if these plates are 
insulators, the temperature gradient must be zero. In general, there- 
fore, it may be expected that six boundary conditions will be satisfied 
Since there are six homogeneous linear differential equations in six 
variables, there is a characteristic -value problem if the boundary con- 
ditions are also homogeneous. Let 




R, 


a, 



6 ( 55 ) 


represent a complete system of six solutions, and let the boundary 
conditions be 
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2a, 


Z, j = 0 


at y = y k , 


k = 1, 


2, 


6 (56) 


where L k is a homogeneous linear function. Then if the solution is 


6 

Z i = ^ A j Z ij(-f) i = 1, 2, , 6 

1 = 1 


there results 




A j Ij k | > 2„,(y k ) 1=0. 


k = 1, 2, ... , 6 


(57) 


(58) 


j 


Hence, there follows the secular equation 


L Jk ( MS , , E, a, c) 

]7 2 


= 0 


(59) 


where 


kjk - ^k j 2 2 j(y k ), * • 


' Z eJ (y k) j > 


t) , k , 2 y , . . , fa 

If equation (59) can he solved, for c, there results 


J 


(60) 
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The condition ^ = 0 gives the limit of stability. 

For incompressible fluids without the effect of gravity, plot the 
curve ci(a,h) =0 in the a-R plane. Here, it has to be done for a 
series of values of M £ and l/F" . 

Conti nuous characte r istic values .- In case one of the conditions 
(56) is absent (ci . the case of "supersonic disturbances" in a boundary 
layer, sec. 5) , no such relation as (6 l) exists, and a solution satisfy- 
ing the remaining five boundary conditions (and certain other conditions 
of boundedness) can always be found. This is the case of "continuous 
characteristic values." The physical significance of such solutions 
will be discussed as the case turns up. 


3. Solution of the System of Differential Equations 
by Method of Successive Approximations 

The exact solution of the system of differential equations (19) to 
(23) or rather (30), (31)? and (33) to (36) is almost impossible. With 
the appearance of the small parameter l/ctR it seems desirable to use 
the method of successive approximations. The general plan of solution 
will be exactly the same as in the case of an incompressible fluid (ref- 
erence 1). Two methods of solution are possible, the first using con- 
vergent series and the second using asymptotic series. 
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ilL-gp j-ution means of convergent series. - In the first method, 
introduce the parameter 


c = (c®)” (64) 

and the new variable 

n “ (y - y c )/e where w(y c ) = c ( 65 ) 

The equations (30), (31 ) , and (33) to ( 36 ) then take the following forms: 


dZ-, 


ar t 


6 


( 66 ) 


dZ. 


dr| 


€ 


(67) 


^■3 _ . „ p‘ ,, . / / M 2 1 

-ieZ^ - _ € 4.3 - . 1 e(w - ci Z 4 + Z = l 

*1 p \ p T 7 


(68) 


cLZg 3. 


dn 


-—7 ’ P j iCw - °) z i + w Z 3 I + iZ 4 | + e 0 (l) ( 69 ) 


i € Hj. 1 - L 
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dZ 6 

dr} 


a f 
( P 



c)Z . + T 1 Z_ 

7 V o 3 ! 

J 


7 - i r , -> -v 

i p * 3 - i(w - c)lfo I { + 0 ( 1 ) 

7 L *J J 


(TO) 


dZ 4 

drj 


- i f | M 2 Z 4 -Zg j- i<x 2 (w - c)Z 3 + a£ 3 0(l) 

1 + ia 2 e ?J _ZT((i 2 + 2(1 i)-7-~-_.. 7 M 2 
J P 


(71) 


The average quantities v , p, T , and So forth, are to be regarded as 
expanded in Taylor’s series in the neighborhood of w = cj thus. 




v 


c = v G ( £Ti ) +_^ (etj ) 2 + . . . , 

2 ’ 


P = P r + pLUi\) + . . . , 


} (72) 


(In general, attach the subscript c to denote quantities at the criti- 
layer where w = c.) The coefficients of the system of equations 
( 60 ) to (71) are therefore convergent power series in e so long as the 
power series (72) are convergent and the condition (37) is not violated. 
An attempt can then also be made to obtain a fundamental system of solu- 
tions as power series of e. A consultation of equations ( 66 ) to (71) 
shows that the solutions should be of the following forms: 
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Dropping t:he su P ers cript zero and recalling the definitions (29) and 
( 73 ) gives, to the initial approximation. 


f X x (p), cp = ex s (n), 0 - X 5 (p), - X 4 (p) 


(80) 


which satisfy the equations 




i w c _ dXx 


Tl 

dr] 3 u lc 


d 2 Xg iv r _ a 


dp 2 Pi c 


dp 


« 

c u c 


F^UicTc 


CX0, 


p Ag = 1 JL C 

Pxc v 


dX, 

7-* = - iA X 
dp 


5 fa c - 7 —l ^2 6 } 


7 Pc ' 


1X3 


(81) 

(62) 

(83) 


where Pj is the kinematic viscosit; 
the Prandtl number. 


coefficient ii x /p, and 


The higher approximations gi vo noriiomogeneous equations, which are 
too complicated to be written out in detail . The homogeneous parts of 
equations are the same as ( 74 ) to ( 79 ) with X'jp' replaced by 
x i (i = 1 , 2 , . . .,6). The inhomogeneous part consists of functions 
of lower orders, and is therefore known. Thus, if the equations ( 74 ) to 
( 79 ) can De solved, these equations for higher appr oxiiua t i one can all 
be solved by means of quadratures. 


If 


i 5 is eliminated from (8l) by means of (82) and (83), a 


dif 1 erential equation of the sixth order is obtained for X 


will give six independent solutions .• 


The corresponding function a 3 


which 
X 2 

and X 5 can then oe obtained from (83) and (82). This is fairly com- 
plicated; fortunately, the case where the Froude number is very large 
is of interest. Thus, as an initial approximation, (8l) may be re- 
duced to 


da 3 


. w 0 dX x 

1 — p — ~ 
t>iC ' dp 


0 


( 84 ) 


the solutions of which are 


= / <H S/ 3 ! 3^] X 12 

/ *- f 

'-J - 


fs^ 

! 1 /3 


D i 

-J 


2 a 




X 


I 3 ’ 


(85) 
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where 


*■( 


„ * \3 / 3 

£cl._ 

Vic 


\ 

\) 


7 


(86) 


Apparently, three solutions are missing. These can he supplied by 

X H ~ J ‘15 = °> X 16 = 0 (87) 

Tne reason these . solutions are , not trivial will he clear when the cor- 
responding 1 unctions X 3 and X 5 are worked out below, 

lrom ( 03 ), it is c.i.ear that the functions corresponding to 

are 


X 


3 1 


/ w c 

- W TT 2 -/ 

\ ic/ 


/id) f 

1/3 j ^ v ^ 


t/^ 


1 (it)t 


t s dt -y Hi *3 |f (U) 1 




u 3 2 


V 1 C ' 


3 *1 1 


/y i \l A5f n / \ f 

iaf ~ ; £ / E ^[ f fHF jp 4 £ 


'K 1 ) 1 2 

3 

h [ I w* pu (f 


(3) 


373 


J 

;V> 


(88) 


x 


33 


- e 


v y ic / 


The functions Zj corresponding to ( 87 ) are 


^ ~ ^3 5 “ = 0 


(89) 


Corresponding to each function X 3 , there are two particular 
integrals X 5 obtainable from ( 82 ) . These are expressible in quadra- 
tures involving Hanlcel functions, for the left-hand side of ( 82 ) has 
the solutions 


X, 


1 

F 


H 


( 1 ) 

1 / 3 


(i£)® 


V°c 


'56 



Indeed, the sets of functions (X 15 , 
two sets of solutions of ( 8 l) to ( 83 ) 
( 87 ) is not trivial. 


•^ 35 * -^ 55 ^ (^- 16 > ^36’ -^ 56 ^ form 

, and it now becomes clear that 
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Next, the asymptotic solutions will he studied, and the convergent 
solutions will he discussed later when the houndary conditions are con- 
sidered. 


(2) §- olut i °n by means o f .asymptotic ser ies.- Analogous to the in- 
compressible case, two asymptotic solutions are obtained hy the most 
naive method of expanding the solutions in powers of (aR)-i. In the 
present case, the initial approximation gives the inviscid equation ( 54 ), 
which is of the second order. The equations for successive higher ap- 
proximate are inhomogeneous; the homogeneous part is of the same form 
as (pi), while the inhomogeneous part is related to functions of lower 
orders. Hence, the integration of all the differential equations can be 
done in terms of quadratures, as soon as (54) is solved. 

Four other asymptotic solutions are obtained by putting 


in ( 30 ) to ( 36 ), where 


Zi 


"i exp f (aF) ; 



( 91 ) 


f i = f'i 0) + fi l) (aR)' , 2 + . . . , i = l, 3 , 4, 
fi » fi 0 ) (aB)^+ f l l} + fi e) (oB)"^ + . . ., 




> (92) 


= 2,6 


} 

J 


while g is independent of (dR) . The initial approximations are 


(Z>. >Z 3 ,Z 4 ,Z 3 ) = (l, 0,0,0) exp f± (aR)' 5 [ /^- (w - J) ay l 

V j 


>(93) 


(Z 2 , Z 6 ) = (aR)'- (± / A- (v - c), OJ exp <±(dR)'~ 


,y 


r n 




~ c)d yj 


(Zi,Z 3 ,Zi ,Z 5 ) = (0,0, 0,1) exp f ±(dR) 


y 


X 


y c 


£ (W - c) dy ; 


V(94) 


(z 2 ,z 6 ) 


- (dR)*- w - c) J exp | ± (dB.) : 


I 


J y< 


. c )dy l 

V 1 Jj 
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ftach of these seta contains two solutions. From these expressions^ it 
appears that the solutions are multiple -valued. Actually, they are valid 
only lor certain regions of the complex plane determined by comparing 
them with the asymptotic expansions of the convergent solutions (85) to 
( 2 °/ • Analogous to the incompressible case, the asymptotic expressions 
hold when (cf. equation (5*^) of reference 1} 

- ~ < arg (£)< g, and - I? < arg (£a c V3 )<g (95) 

simultaneously. If e is very close to a real number, this means that 
the expressions (93), (94) represent solutions in a connected region 
which contains at least a substantial portion of the real axis. This 

fact will be seen to be of significance in discussing the boundary- value 
problems . 

Similar considerations hold for the solutions of the inviscid equa- 
tions . . These solutions, appear to . possess a logarithmic singularity at 
he point w = c. As in the incompressible case, the asymptotic expan- 
sions of the convergent solutions bring these solutions into correspond- 
ence with X 33 ana X^ , and the restrictions (95) explain this appar- 
ently multiple-valued nature of the solutions. 

. A f aiogous to incompressible case, there are points on the real 
axis where the asymptotic solutions fail in the case of damped dis- 
turbances. These are interpreted as "inner viscous layers," where the 
effect of viscosity is not negligible no matter how large the Reynolds 
numoer is. Considering the conditions (95), it is seen that there are, 
m general, f our of them. For the lines 


arg (£) « - arg (fj = | 


arg (to c l/z ) = - ll t a rg (£a c l/5 ) c £ 

0 6 

intersect ^the real axis in four points, if a is not a constant and 
c i < 0 . These are the points where the inviscid solutions fail. The 
four points reduce to the single point y c when c is real, and there 
is no intersection when c± > 0. Hence, the four inner viscous layers 
coalesce into one in the- case of neutral disturbances and disappear 
completely for self-excited disturbances. The significance of these re- 
sults for the studies in part II is also similar to that in the incom- 
pressible case. 


33 



BACA TU Bo . 1115 


4. Boundary-Value Problems 

Having obtained the solutions in convergent series of € and in 
asymptotic series, the boundary-value problems discussed briefly in sec- 
tion 3 will now be enlarged upon. The physical requirements give rise to 
certain mathematical conditions on the real axis of the complex y -plane. 
In general, the boundary conditions for the velocity disturbances are 
independent of the temperature disturbances. For example, for flow in a 
channel with walls at yi and y 3 , the boundary conditions 

ffai) c <P(yi) = f(y 3 ) * <p(y 3 )--= o ( 96 ) 


must be Satisfied whatever the conditions on the temperature disturb- 
ances may be. These conditions are identically satisfied by the solu- 
tions Xi 5 and Xfg , to the proper degree of approximation. Thus, the 
quantities defined by ( 60 ) vanish if j = 5,6 and k = 1, 2,3,4 

(say). Hence, condition (59) reduces to 



1 , 2 , 3, 4 


(97) 


-lj^^ggg_fcgyi 8tic -value problem therefore does not explicitly depend on 
the, tempe rature disturbances in the initial approximation . Indeed, after 
the temperature disturbance corresponding to the characteristic oscilla- 
tions has been determined, it is always possible to satisfy the boundary- 
conditions for the temperature disturbances by including a* suitable lin- 
ear combination of the solutions X 55 and X 56 . The corresponding ve- 
locity disturbances are identically zero and will therefore not interfere 
with the boundary conditions imposed upon X 55 and X 56 . Such a situa- 
tion is in general the case. The characteristic-value problem therefore 
becomes very similar to that in the case of the incompressible fluid. 

Two inviscid solutions f li? , <Pj.p are derived from (54) and two viscous 
solutions from ( 85 ) and (88): 




r* 

2 

ij 


s 



^( 98 ) 



As explained in the last section, the asymptotic solutions hold in a 
connected region containing most points of the real axis except a 
neighborhood of the point y c . Thus, there will be no difficulty in 
using the solutions fj^,-., 9 i ,2 for the discussion of boundary - value 
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problems . The path leading from one boundary point to another lies es- 
sentially in the lower half plane. 

Now, let these results he applied to the case of the boundary layer. 
Begin by making a careful investigation of the boundary conditions. At 
the wall the conditions f = cp = 0 hold. However, f my be replaced 
by a linear combination of cp and cp* . In the case of the viscous so- 
lutions, 93,4 + i^ 3,4 B 0, so that f 3>4 = icp 3 4 . For the inviscid so- 
lutions, -the relation * ’ 


f 




T - M'-(w - c ) 2 


(99) 


holds (equations (49) to (53)). At the wall, w = 0, and 


(99a) 


Thus, the condition that f vanishes at the wall may be replaced by the 
condition that a linear combination of <pi, 2 /p 3 4 , and <p x ,g vanishes. 

Analogous to the incompressible case (cf. equation ( 6 . 8 ) of refer- 
ence x) , the condition of boundedness at infinity rules out the solutions 
x 4 and cp 4 . it is convenient to take the lower limit of integration in 
f 3 and CP 3 a t + c °. Then, these solutions vanish rapidly as y — ~ s»oo. 
Thus, for y » 1 , the inviscid solutions dominate. In the incompress- 
ible case, the inviscid solutions behave like e ±ay . The condition of 

boundedness therefore leads to Cp ~ e This is conveniently express- 
ible as 

cp* 

~ + a — > 0 as y — > o» ( 100 ) 


In the present case, a corresponding condition must be established. 
However, the situation turns out to be more coniDlicated. Consider the 
equation 


J± f (■¥_ - c) CP* - w*cp 1 _ q 2 (w - 
\ T - M 2 (w - c) 2 J T 

obtained from (54) by dropping the term in \ 

p- 

more in detail in part II, the behavior of the 
is given by 



( 101 ) 


As will be discussed 
solutions as y — $» a, 
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9 - e"^ 7 , if 3 = a Jq = a/l - W (l - c ) 2 4 0 


( 102 ) 


where j3 is uniquely determined if a cut is drawn along the negative 
real axis of the complex Q-plane. Then it is clear that the real part 
of 3 is always positive if Q is not on the cut, and hence the solu- 
tion 9 ~ e + ^ y must be rejected. Thus, there results the condition 

.3L + 3 — > 0 as y — « ( 103 ) 

This reduces to the incompressible case if M = 0. If Q lies on the 

cut (c = C r < 1 - g j, the situation is more complicated. The situation 

will become clearer only after a thorough investigation of the inviscid 
problem (sec. 7 ). 

Except in the case 1 - M'(l - c) ^ 0 , the characteristic-value 
problem is therefore almost the same as that in the incompressible fluid. 
The characteristic values are given by the determinantal relation 


where cp^ 


9xi 

921 9 31 

* 2 ^ 

t 0 » 

Tx9xx + M cq/xx 

Tl92X + M W) CCp21 * 

T , ..“* **■ /*rv 

T x - M a c 3 

T, - M 2 c 2 31 


‘ „ * 

912 + j3cp 12 Cp 22 + |Bcp 22 0 

and cp 2 are any two linear independent solutions of ( 101 ), and 

<Pij = <Pi(yj), <pl j = 9 i(;'j), 1 ,J = 1,2 (10k) 


y 2 being the coordinate of the "edge" of the boundary layer. Strictly 
speaking, the value should be substituted for y 2 . However, for 

J > Jz> the solution of ( 101 ) is practically identical with e ± ^^ r . 
Tims, it is a good approximation to impose equation ( 103 ) for y = y 2 . 
Naturally, the larger the thickness of the boundary layer is taken, the 
better is the approximation. 


The determinantal equation may be written in the form 


E(a,c,M a ) -= E(z) 

36 
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where F(z) is the function of Tietjens (reference ll) 


F(z) 

with 


-z 



+CO 


H 


(i) 

1 / 3 



3 







(106) 


(107) 


rjx being the value of ri at y = y x , the solid boundary. The function 
E(a,c,M c ') depends only on the inviscid solutions, and is given by 


(yi 


y c )E(a,c,M 2 ) 



* 


Ti9ii + M 2 wic 9 n , 


9ll 9X2 + 39X2 


T x - M s c 2 912 + 13912 


» 

*7* 

T i92i + M 2 w x c9 2 x t 


921 922 * 3922 

1 

Tx M'V 2 922 + P922 1 


( 108 ) 


The manner in which the viscosity coefficient enters the final equa- 
tion (106) is noteworthy. As compared with the incompressible case, it 
amounts only to a change of the definition of z. By referring to (106), 
(6k), and (6p), it is seen that this amounts to the replacement of R by 
^/ u iC‘ This means that the Reynolds number defined in terms of the free 
stream v elocity and the kinematic viscosity coefficient at the critical 
layer (instead of that in the free stream) is the quantity governing sta- 
bility phenomena. This point must be kept in mind whenever it is neces- 
sary to compare a case of homogeneous temperature with a case of inhomo- 
geneous temperature. Greater detail will be given in discussing the 
stability problem in a real fluid (pt. III). 

In the case 1 - M s (l - c) s < 0, it is possible, of course, also to 
impose the boundary condition ( 103 ) , with 0 imaginary . The same equa- 
tion (108) holds. But the general discussion of the physical signifi- 
cance of the solutions is more complicated. It will be discussed more in 
detail in part II. 


The inviscid case . In the limit of infinite Reynolds number, the 
the relation (106) becomes 


$11 ‘Pi a + 09i2 


921 9,22 + 09? £ 
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= 0 


E (a, c,!! 2 ) = 0, 


or 
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Thle corresponds to a solution of (101) with the boundary conditions 

cp(yi) - 0 , <p* (y 2 .) + (3cp(y 2 ) « 0 (HO) 

Consideration of this boundary problem gives the asymptotic behavior of 
the relation ( 6l ) in the form 

c = c(a,M s ) (m) 

This will be discussed fully in the next part. 
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II - STABILITY IN AN INVISCID FLUID 
5. General Considerations 

It has been shown that, in the limit of infinite Reynolds numbers, 
the problem can be treated with viscosity neglected, provided proper care 
be given to the inviscid solutions. Such investigations will naturally 
give some information to the stability problem in a viscous and conductive 
fluid. Indeed, the complete calculation of characteristic values, in par- 
ticular, of the limit of stability, can be carried out, once the invisci^ 
solutions are known. It is therefore advantageous to study the inviscid 
case as a prelude to the actual case, with the expectation that certain 
important characteristics may be obtained. 

However, it must be noted that the results obtained in this case can- 
not be applied directly to the viscous case without modification. Thus, if 
only stable (damped and neutral) disturbances can exist for a given flow in 
3X1 ' inviscid fluid, it cannot be concluded that unstable disturbances cannot 
exist under the action of viscosity. However, if unstable disturbances ex- 
ist in the inviscid case, the flow will still be unstable when viscosity is 
taken into account. For if the continuous dependence of c^ on R is con- 
sidered, it is evident that c^ cannot remain less than or equal to zero 

for all ignite values of R and still become positive as R becomes infi- 
nite. 

This investigation will begin with a careful study of the analytical 
nature of the solutions, especially for y becoming infinite. It is 
found that the disturbance there takes the form of progressive waves out- 
side the boundary layer. For more detailed discussions of their properties , 
it is found convenient to classify the disturbances as "subsonic, " "sonic," 
or "supersonic" when the x-component of the phase velocity of the disturb- 
ance relative to the free -stream velocity is less than, equal to, or 
greater than the mean speed of sound in the free stream. 

The amplitudes of these waves go to zero as an exponential function of 
the distance from the solid boundary, except in the case of neutral super- 
sonic disturbances. To an observer moving with the velocity of the free 
stream, the waves are propagating opposite to his direction of motion for 
neutral subsonic disturbances. For a general disturbance, the direction 
of propagation is inclined outward if the wave is amplified and inward if 
it is damped. For the neutral supersonic disturbances, there may exist 
both an incident wave and a reflected wave with (in general) non- vanishing 
amplitudes at infinity. 

Analogous to the incompressible case, an attempt is made to establish 
necessary and sufficient conditions for the existence of certain types of 
disturbance. But a consideration of energy relations is found to be ex- 
tremely helpful. This is carried out in section 8. 

In the case of neutral or slightly non- neutral subsonic disturbances, 
the physical situation for the compressible fluid must be quite similar to 
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the situation in the limiting case of an incompressible fluid. There- 
foie, it should be possible to obtain a general criterion for the exist- 
ence of slightly amplified subsonic inviscid disturbances, analogous to 
the Rayleigh- To llmien criterion for an incompressible fluid (sec. 9b). 
After such a criterion is developed, mean velocity- temperature profiles 
could be readily classified according to their relative, stability at very 
large Reynolds numbers, and the effects of the compressibility w.nd con- 
ductivity of a gas on the stability of laminar boundary layer flow can 
be evaluated (sec. 11). 

In the case of neutral supersonic disturbances, both incoming and 
outgoing waves may exist, with the amplitudes of the incident and re- 
flected waves unequal in general. Except in the particular case of a 
pure neutral outgoing or incoming wave, there is therefore no character- 
istj.c- value problem^ or rather, the characteristic values are continuous, 
and not discrete. By utilising the results of the investigation of the 
energy balance for a neutral inviscid disturbance (sec. 8) a general ex- 
pression will be obtained for, the ratio of the energy carried out of the 
boundary layer by the reflected wave to the energy brought into the 
boundary layer by the incident wave (sec. 10). With the aid of this ex- 
pression for the "reflectivity," at least a necessary condition for the 
existence oi a pure neutral outgoing or incoming wave can be determined 
(sec. 10). 


6. The Equation of Inviscid Disturbance and the Analytical Nature of 

the Inviscid Solutions 

-k- 1 ^4e limiting case of infinite Reynolds number and infinite 
iroude number, the disturbance equation for <p reduces to the follow- 
ing linear differential equation of the second order (cf . (101)). : 


d f (w-cjcp'-v'cp 
dy V T-M a (w-c) 2 


tp 


or, in the self-adjoint form, 


a 

dy 




where 


(y) -r 4 ( y ) „ J_ L(t 

J V-C dy dy 


( 112 ) 


( 113 ) 


(114) 
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There is also the relation 


* = .. lyM 2 f (v--c)9 t -w , 9 ~1 = iv q> r + if 
P l T-M 2 (w-c>J ■ v-c 

The second part of the last equation may also be written as 


Tcp» -MSw'(w-<;)q) 
T~M' (w-c) 2 


( 115 ) 


(116) 


omce the coefficients of the differential equation (l!2) are entire 
functions of the parameter a 2 , the- two particular integrals cp x and' <p 2 

f£ n J? ls ® 9 ^ atl ° n J must &ls ° be entirs Unctions of a? . Series develop- 
menuo ol cp x and cp 2 in powers of a 2 are therefore uniformly con- 
vergent for any finite region of a, for a fixed value of y. except 
when y is a singular point of the differential equation.* 


If the series development cp = ct/° ^ + a, 2 cp ^ 1 ^ 
substituted into (l!2), two particular integrals ©■» and rpo 
by successive quadratures. 


. . is 

are obtained 


00 


<Pi (y; a 2 , c, M 2 ) = (w^c) _) a aJ % 2K (y, c, M 2 ) 


n=o 


^2 (y; eh } c, M 2 ) _ (w— c) a^^gn+i (y, c , M 2 ) 

n=o 


( 117 ) 

(118) 


where 


fcan(y; c, vt) - / 


," y --M 2 ) dy r 

J 7x ' (& - 0,E " --'jx T 


1-0 (yj C, M2) = 1 


x l 2 n- 2 (y; c, M 2 } dy 


f ( 119 ) 


It will be shown later (sec. 7b) that the point- at which T equals 
M 2 (w-c ) is only an apparent singularity of (112). The point y = y 
and the point at infinity are the only true singularities. 
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and 


k2n+i(y;c,M'') 


rJ 


r 


1 n 7 (v-g)' 


M } dy 


L L(w -° )2 J T 


k i (y; c, M 2 ) » P f— - 


- M a l cb 


, / '^(w~c ) 2 J 

'yi 


x (y;c,M 2 )dy, n> 1 


\ ( 120 ) 


In these integrals, the lower limit is taken at the wall (y = yi ) merely 
for convenience. 

In order that cp^ and cpg may he valid approximations to the regular 
solutions of the complete disturbance equations (.19) to (23) all along ’the 
path of integration between the points y = y x sna y = y on the real 
axis, that path must lie wholly in a region in which the asymptotic ex- 
pansions of the regular solutions for large values of cffi are valid. 

The asymptotic expansions of these solutions for large values of <xR 


hold only in the range defined by 



< 


arg. 





(cf. (95)). 


Consequently, the path of integration between y x and y must be taken 
below 1 the point y = y c . (See fig. 1.) 


Contrary to the statement made by Tollmien (reference 4), the proper 
path must be taken be low the point y = y c regardless of whether this 
point lies above, on, or below the real axis. This question, which has 
been responsible for a certain amount of confusion, was finally clarified 
recently by C. C. Lin (reference l). 
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It is now possible to define a region in the complex y-plane in which 
the solutions cp x and cp 2 are everywhere analytic in the variable y 
and the parameters cP , c, and M 2 . Consider the simply connected region 
t' ^ nc lhdes the end-points y = y 3 and y = y 2 but not the (singu 

point y = y c (fig. 2). The region R* and the region S’ in the 
neighborhood of the point y = y c can be made mutually exclusive. Pro- 
vided w(y ) f 0 in the range y x < y < y 2 , the relation c = w(y) maps 
the regions R* and S’ in the y-plane into the mutually exclusive re- 
gions R and S in the complex c -plane. 


to 

y 


R' 


the coefficients of 


and the parameters 


co 


If y is now restricted to R* and c 
(112) are analytic functions of the variable 

and bP, and the solutions cpi and cp 2 must also be analytic func- 
tions of (yi, a 2 , c, M* ) . So far as the characteristic-value problem 
is concerned, the 'analyticity 1 of the solutions cp x and 9 2 in a simply- 
connected region enclosing the boundary points y 3 . and y 2 is assured*. 
Unfortunately, this argument fails when c = 0, 
point y = y c , w = 0 coincides with the point 

boundary, and the regions R* and S' cannot possibly be mutually ex- 
clusive. This special case will be discussed briefly in sections 9a and 
9b. 


because the singular 
y = yi at the solid 


7. Further Discussions of the Analytical Nature of 


the Solutions; 


Their Behavior around the Singular Points of the Differential Equation 


Although the analytical character of the solutions <p x and <p 2 in 

+V* 3 v G f i0 ? R , ls of great importance for the characteristic-value problem, 
l®l ehavior °;f and ep 2 in the neighborhood, of the singularities of 

1 • t S , lffi P° r tant in the investigation of the physical mechanism 

oi instability. 


( a ) Singularity at the Point 


w 


c: 


Uhless the quantity. 


d 
Lay 


( 


j 


w=c 


I a (v*\ 

Ldy VfV 


w=c 


vanishes. 


the point y in the complex y— plane, is a regular singularity of the 

differential equation (112). Since (w-c) and t' are analytic functions 
y everywhere in the finite region of the complex y-plane, they can be 

Pn In ' I 'Cl tt ! 4 _ nil . t i ' K 


of 


~ — > .i. UJ.IC/ 

developed in Taylor’s series around the point 


y 


y c = f; )i as follows: 
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V C .• . s> 


w_c = w c * (y~y c ) + (y~y c ) + 

21 


( 121 ) 


T c " 


T T c + T c* (j-y c ) + — — (y~y c ) + . . . 


( 122 ) 


Upon substituting the series developments (l2l) and (l?c>) i n + n 

sss g ' iyrw 4 


two 

the 


9i = (y-y c ) g(y-y c ) 


(123) 


T„ 


9a = g2(y-y-) + 




a (v) 


<37 


W=C 


91 log (y_y ) ( 12 k) 


where Si and g 2 are analytic functions of (y - y _) „ 2 c v 2 

and g, (0) = 1 -la «. /7\\ m / , / * ^ ana & , 

he suhlected +-,"+>» * & ( ) - T c /w c ^ o. These two solutions must 

(113) J SL SaJ2e restrictions (?5) as the solutions (117) and 

i n * Con8e< l uen tly, in passing from RZ (y - y r ) > 0 to KZ (v v ) 
< 0 the correct path lies below the point y 1 ?nd tl- niSnZ/® 

analytical continuation of (124) for (y - yj < ^°% 8 1 P P ® r 


92 = g 2 (y-y ) + ,r d i /V,N ) 1 
U J C ; - J 

c c 


9i 


log | 


V— 


- 1 


ill (l2ha) 


solutions Physical problem, of course, only the properties of the 
solutions and along the real axis are important. If 0 > 0 

“sturtenoe) the point y = y 0 llee above the real J ^ 

the- solutions are regu lar along toe real axis. In this case, the effect 

/ l3 l In ° e f is related to 9 by (116), the discontimtotyluffered by”" 
In from (y-y c ) < 0 to (y-y c ) > 0 leads to a phase 

trSS^oJtne^frorJb 1 * U ^ phaee shift " hloh Possible «*> 

(sec. 8). 87 f the mean flow to the disturbance, or vice versa 
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of viscosity and conductivity on the disturbance is negligible in the 
interior of the fluid for very large Reynolds numbers. However, if 
c i < 0 , the inviscid solution (124) cannot possibly be valid all along 

the real axis (fig. l). If c^ = 0 (neutral disturbance), there is a 

critical layer of fluid at the point w = c in which the velocity varies 
very rapidly (f log { y_y c f ), and in which, therefore, the viscous 

forces must be taken Into account even when the Reynolds number becomes 
indefinitely large. If c^ < 0, (damped disturbance), there are four 


such Inner critical layers, because the lines 
(fig. l) and the lines arg w h <y-yc)} 



7n n 

-, which delimit 

6 6 


the region of validity of the solutions <p x and <Ps, intersect the real 
axis in four points. 


With the aid of the equations of motion and the relations ( 115 ) and 
(ll 6 ), the physical situation in the neighborhood of the point w = c 
can be made still clearer. It is not difficult to show that the rate of 

change of the quantity 1 p*(* where £ * is the vorticity, for any two- 
dimensional motion in an Inviscid, non-conductive compressible fluid is 
given by the relation: 


JL (p*t*) = H* + ~ d( p **P*> (125) 

dt* ~ &t* p*' y*) 


In the present case, if c = c r , then from (125) and the equations (49) 
to (53), there is obtained 

ipr (p*f*) = 0 at w = c ( 126 ) 

or 

<P — (pw*) =0 at y = y c (127) 

that is, the transport of the quantity pw ? across the plane w = c 
must vanish. It will be shown later (sec. 8 ) that it is impossible for 
^(yc) = 9 C to vanish if cp(y) is a solution of the disturbance equation 

1 * 

The quantity p*£* is related to the density of angular momentum. 
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(112) which satisfies the boundary condition cp = 0 at the wall. If the 


value of 


c = e. 


is chosen so that 


L dy ^ PW ' • 1 t 
t ^ -j w=c 


then the trans- 


p ort of pw across the plane w - c can be balanced only by the diffu- 
through the action of viscosity. It can therefore be con- 
cluded that a neutral disturbance free from the effects of viscosity in 
he interior of the fluid can exist only for velocity-temperature profiles 


for which 


_d_ 

dy 


(pw*) 


_d_ 

dy 


= 0, at some point. 


From the energy equation ( 52 ), 
(124a) , it appears that if c = c r 
comes indefinitely large as w -■> c , 


the relation (ll6), and (123), (124), 
(neutral disturbance), then 6 be- 
y->- y c . Even if the quantity 


d ^w*\ 
dy ) 


JC 


vanishes, the conductivity of the fluid cannot be neglected 


in the vicinity of the point w = c unless T‘ = 0, which is not gener- 
ally the case. However, the mathematical results obtained in uart I in- 
dicate that the influence of the conductivity on the "viscous"’' solutions 
oi the velocity components is only secondary for Keynolds numbers of the 
order of magnitude of those encountered in most aerodynamic problems. 

(*>) Singularity at th e P oint T =. M s (w - c) g ; 


In the case of the neutral supersonic disturbance, outside the 
oundary layer, the relative velocity between the mean flow and the x~ 
component of the phase velocity of the disturbance is always greater than 
the mean sonic velocity. At some point within the boundary layer, the 
relative velocity must be equal to the local mean sonic velocity a. 

v 2 m 


Since 


(%,. 

Vtrv 


~~ this point will be reached when (w— c) = 


a* 


u 


T 


M 


Although l (y) — > 00 as T — M 2 (w-c) 2 , by means of a change in 


dependent variables it is not difficult to show that this point is only 
811 singularity of the differential equation (112). If the de- 
pendent variables are chosen as cp and -, rather than cp and cp‘, 

then by utilizing (112) and (ll6), a new system of linear differential 
equations of the first order is obtained: 
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dcp i . T - M (w - c)^ jr w* 
+ 

dy . (v - c ) p w ~ c 



dy 


IrM'V 


(w _ c) 




(128) 


(129) 


The only" singularities of equations (128) and (129) occur, at the point 
w = c and the point at infinity. So 'far as the disturbance is concerned, 

the physical significance of the point T = M 2 (w ~ c) 2 lies only in the 
fact that it marks the point of transition between the supersonic and sub- 
sonic fields of flow. 


For. 

the neutral sonic disturbance ( 1 — c 

V 

M 

the 

point at which 

2 

c) 

, T 

(T 



equals — moves out to infinity; 

1, w - 

1, as y -> co ) . 


M* 


The physical and mathematical problem is more difficult to investigate 
in this case, because equation (l!,2) has an essential singularity at in- 
finity and the asymptotic behavior of w(y) and T(y) as y -> cc ± B 
somewhat complicated. In the next section, the behavior of the sonic dis- 
turbance as y — > co will be discussed in some detail. 

(c) Behavior of the Inviscid Di s turba nce as y — v 

Boun dary Conditions and the Charac te ristic-Value P rob lem 

As y —> * , T — 5» 1, w — •$> 1, w*-> 0. If a 2 (l — M 2 (l~c) 2 l f 0 
the disturbance equation (112) takes the limiting form 


q> = a 


/ 


M a (1 - 0) 


2 


J3y 


Equation (130) has the solutions e ’ " and e 

*2 


J 

+‘3y 


<P 


(130) 


, where 8 = a ,/ £1 

and 0=1- M‘"(l - c) 2 . It follows that the equation (112) has a funda- 

:t{3y 

mental system of solutions behaving like e “ as y » . f 0 define 
8 uniquely, it is necessary to introduce a "cut" along the negative real 

axis of the complex Q —plane. Begardless of whether > 1 _ M 2 (l-c ) 2 '1 - 0 

^ x r * ' /P > 

Ki 


J 
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(subsonic, sonic, or supersonic disturbance) 1 the real part of* p will be 
so long as - « < arg (,q)< it. Since the physical conditions of 
the problem require that cp must be bounded as y <* , the solution 

+Py 

e must be rejected. Therefore, this solution cp must behave like 
— Py 

e as y -> oo . 


Solutions of the type e , when combined with the factor ' e i«(at-«t) 
evidently represent progressive waves, but it is necessary to be careful 
in discussing its direction of propagation. A disturbance which is pro- 
pagated outward with respect to a fired observer at the wall is actually 
an wave relative to an observer moving with the velocity of the 

mean flow outside the boundary, layer, and vice versa. This fact can be 
readily appreciated by referring to. figure 4. The wave fronts moving 
ogfcff gft and- also downstream with a velocity c r relative to the wall are 
overtaken by the observer moving downstream with the velocity 1 relative 
to the tall. To such a moving observer, these wave fronts appear to be 
propagating inwar d < -the situation is obviously reversed. -for the wave 
fronts moving inward and downstream with respect to the wall. The same 
conclusion can be reached, of course, by referring to the analytical form 
of the disturbance. A wave front moving outward with respect to the fixed 

wall will have the form However, for the observer 

moving with the free stream velocity, the x-coordinate is x* = x — t, 

and the wave front has the form e ia ^* + ^ '” c ' ^©^ i n his (x 1 , y) 
coordinate system. The wave front is obviously moving inward in this* 
system. If > 0 (disturbance increasing with time )7 Then Q 4 > 0 

and p i > 0; the disturbance takes the form of an out going wave of ex- 
ponentially damped amplitude (in y) as y -> a> , if c < o (disturb- 
ance damped with time), then Cl ; r < 0 and p^ < 0; the disturbance takes 
the form of an incoming wave of exponentially damped amplitude (in y) 
as y-^ oo. if Cl = 0 and Cl > 0 (neutral subsonic disturbance )“ 

the disturbance is propagated parallel to the x-axis, and the amplitude 
is exponentially damped in y as y — > a> . Thus, for — it < arg (q)< it, 
the boundary condition at y = y s is cp f (y) + pcp(y) = 0, and the* 
characteristic values are discrete (sec. 4). 


l 


The curve 



corresponding to the condition 


1— M 2 (l-c r ) S = 0 


for a 
sonic 


sonic disturbance divides the complex Q-plane into a region of sub— 
disturbances and a region of supersonic disturbances (fig. 3). 


48 



MCA TW Ho. Ill' 


If c± = 0 and Q < 0 (neutral supersonic disturbance), then J~Q 

is purely imaginary, and both solutions of equation (130) ( e ±a v /_ ^ 1 7 A 

f inite^equation^S^l ) ^ °°^ respondi fS pressure disturbances are also 

exist but in wnl'ni ii 111 ’kkis case, 13001 incoming and outgoing waves 
£ a -v 8 ! u 1 they are not of equal amplitude. This phenomenon 

at theTOl! f discrete. A solution satisfying the boundary condition 

?a C f ( fT™ hi7) ^T?,?M f0Und . f0P arMtrar > r values of c and a. In 
* . m (117) and (ll8), such a solution for y i < y < y 2 is 


CO 

<p(y) - (» - o) £ «"k SBHl Cjr» c, M s ) (131) 

n~o 

- 0 ^fneutral ti0n , a i'’ ~ M n ” c ^' } f rj breaks down when l-lf(l_c) 2 
7 C “ 5 “ tr *l 80nl0 Haturbance) or when a - 0. In the latter case, with 

and fnA} 1,1 ~t n0t difficult to verify that the solutions (117) 

the point at infinitv°of a as a 0 even as y -> ® , although 

P infinity of the y -plane is an Irregular singularity of the 

e quation (llg). Indeed, these solutions behave like e~^ y (l 


JT~ M s (l - c ) 
a (1 - c ) 
ing case a -$> 0. 


c) and 


- ( 


-3y +3y\ 

e -e ), respectively, as y > 


in the limit- 


For the case of the neutral sonic disturbance, ( 1 .. M 2 (l - c ) 2 1 = 0 
and the asymptotic behavior of the inviscid solution^ as v ~x> .«*. * 

The asymptotic behavior of the mean velocity w(y) for th^ 

Blaslus SofS^rS? 6 fl ^ d i °' m,iarJ larer 18 SImllar ^t of thf 
xasius profile (reference 11, equation 10, p. 228 ). 


1 - w 


,a ° _Xz 
e dz, 


for 


J » 1> (X « const.) 


( 132 ) 
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the , prandtl number is unity ana the mean free- 
ze ^(y is a Sloue?u2™tf flOT JV®™. «*> »» W,r*l 

J unique quadratic function of w(y). Thus, (reference ll) 


T = T x + J l 1 M 2 


(T x - 1) | w — LjI 

J 2 


,,2 a 

M w 


Tu that case , sine© c = 


c Q = 1 - ± , 
M 


(133) 


T - M^w - c) 2 „ (l - w ) F for y » l 


(134) 


yhere F is a positive constant, 
tj&ke the limiting form? 


The differential equation (112) crust 


d j~ (w - c) 2 , -| a 2 (w - c)‘ 

dy \ (T~- '"w)F ^ j ~ t 




a 


i 0 


(135) 


or, 


TT» 2 

•j/ n , v . # a ' 

^ + = — (1 ~ w) F \|/ 

1 - v T 


( 136 ) 


where t = <P/(v c). If the physical condition that \j/ be bounded as 

" iK£P ° 8ed ' then y ^ „ (1 - v)* = 0, and equation (l 3 6) 
admits of two possible solutions: 


ll/ constant, cp - w - c constant, as 


7 -s» oo 


(137) 


or 


* ' ~ 1 ~ <P(y) ~ / (l - w)dy 0 


as y — > » (138) 
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The asymptotic behavior <p ~ (w- c) implies that 


&y V P 


f It 'N v 2 It 2 ( 

\p / a > 02 ? - ~ a (y + A) as y -$> oo , by equation ( 129 ) ( 339 ) 


'-jo 

<P(y) ~ / (l — w)dy, 

* rr 


then 



0 ) — * — 5> constant 
P 


as y — $» .» 


(140) 


remain^ fteite Lr ? d } 8turl)aace i . , if ^he press-ore disturbance is to 
L infinity (equation (l4o)), then <p must approach zero ■ 

Zst iTrtlrZ* becomes infinite (equation (138)). The solution (137) 
must be rejected, and tne characteristic -value problem may be eroected tn 

Sticrs t ^ 81 ' acts n^4 0 raluss - aa con41 “°-- 

1 ition in this case will be discussed in section 10c. 

but ° nl7 a finite Sradient of pressure disturbance is required 

tions (iSfJd 6 r ^®f sur \ di f^f lce itself to become infinite, both solu- 
continuoul ®. included , a f ld the characteristic values become 

•+ - physical significance (if any) of this solution is not clear 

a ^reasiMo^uM f °”f? Et . ana4 ° SOUB to the caBe of the steady flow of ’ 
looaf ^ « vioinilff of the speed of sound, „here snail 

ve4 cross-sectional area hounded hy stream lines produce 

-y t^ge local changes in the velocity and pressure. 

It should be noted that if a = 0 and c = 1 - l, both solutions 

M' 

(wntwf “ 4 ”° incl i d9d - ( Sm of sec. 9 a.) From the rela- 
as y In tMs b case en ^ preSSUre dist ^ance remains finite 


8. Energy Relations for an Inviscid Disturbance 
T J e disturbances have just been classified into (l) self -excited 

parallel to the 
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toth inward and outward (Q < 0). It is very Interesting to consider 
the energy relations in all these cases. 1 

In the first case, there is no doubt that energy must pass from the 
mean flow into the disturbance, because the amplitude of the disturbance 
is being Increased and energy is being carried to infinity by the wave at 
the same time. In the second case, the opposite is true. In the third 
case, there is apparently no transfer of energy between the mean flow and 
the disturbance. In the fourth case, energy is being carried in and out 
by the waves; whether energy will pass from the mean flow to the disturb- 
ance, or vice versa, depends upon whether the amplitude of the outgoing 
wave is greater or less than the amplitude of the incoming wave. 

For the two cases of neutral disturbances, (3) and (if-), it is possi- 
ble to clarify the physical situation by considering the time average 
over a period (which is well defined for neutral oscillations). 

Since viscosity and conductivity are disregarded, and the neutral 
disturbance Is harmonic both in x* and t*, the average time rate of 
change of the total energy per unit volume over one period and one wave- 
length must be zero; that is. 


dt* at* 


(Ui*) ^ 

2 J 


+ P* __ ( Cp T*) = 0 


(l 4 l) 


or 


2 

d f (u.j*) •'j dp* 

dt* l 2 J dt* ( l2f2 ) 


By neglecting triple and quadruple correlations, and utilizing the dynamic 
equations like (49) to ( 53 ) to carry out certain reductions, the energy 
balance for the disturbance is obtained in the following form: 


p* ■ — s 

dt l 


f ( v>‘ 


— - A* 


u. 


u 


- c)u. * 
j ctoCj* 


u- 


4ft 


dp** 


= 0 


( 143 ) 


These investigations will also form the basis for the discussion of 
the necessary and sufficient conditions for the existence of a disturbance 
(secs. 9 and 10). 
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Btu *’ a f „ 

Now the quantity p*» is equal to - El I jr(q>* + iccf) L where 

dXi ' 2 J 

the symbol ~ denotes the complex conjugate. From the relation (115), 

there results 5. = ±y 9. L±JJL f BO that the quantity in the bracket of 
p V — c 

the above expression is purely imaginary when c is real. Hence, 
dui* 1 

p - = 0. CJsing this relation, ( lh-3 ) can be reduced to 

C£ji 



CtU 

— p Uj * f Ui * 1 

1 J ax.* 

3 


3 

r (%"■’ p*‘) = 0 (i44) 

3 V 


The relation ( 144 ) holds for every point in the fluid. Hence, 



- p* Uj_** Uj** 


clu . * -\ 

1 I 

3x“*/ 


dV = 


M. 


/"Wj 

' / / 




f P** 


H*'/ dV 


pr> 

= // p* 1 %*’ as (145) 

■Jj 


where V is a given volume of fluid, S is the bounding surface, and 
u n *’ is the component of the velocity perturbations normal to S. Let 

V be a rectangular parallelepiped of unit dimensions in the x* and z* 
directions, extending from the solid boundary (y = y x ) "to inf ini tv" in 
the y* direction. Then with the condition v** =0 for y = y x , (ih-5) 
becomes 



_ . — du* 

- p* u* f v* f dy* 

dy* 



du* 

dy* 

dy* 


lim 

= p*» v ** 

7 —> 


(l46) 
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T rgy P 7”>f te J OTtBard , ty the disturbance in unit time 
/ = constant - (j large) Is equal to the total energy 
transferred in unit time from the mean flow to the disturbance by the 

action of the shear- stress t * = — p* i^'v**' within the boundary layer. 

It can be verified that 


| XI (f$) 


(l J +7) 


where cp denotes the complex conjugate of c p. By making use of ( 176 ) 
follo3ng B f°™ uta C0rrela * l0n can 18 bought *»*> the 


2 T ~ M 2 (w ~ c) ; 


Im (cp*cp ) 


(148) 


or rather 




a T 

~ Z~~r:s~, 7? ( < Pr < Pi t - ^r') 

2 T _ M (w — c) 2 - 1 


(149) 


Wien c -- c r , the coefficients of the differential equation (.112) 
are real. If cp satisfies (112), then q= r and cp. must also satisfy 
the equation independently, and the expression in brackets in (]4o) is 
the SiS °s th ° tTO BOlutlons > ^ definition. For equation '( 112 ), 


nt' 

- / - dy r 

W (<p r , cpa ) = e / i =k < T — f 


a , 0 h 

M (w - c)' > 
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from the relation (ll4), where 3c is a real constant. Hence, from (1.49), 




a 


kT, t » 


a 


W' TT* 2 


k 


(151) 


o u o 


Thus, if c = c r , the shear stress is c onstant wherever ¥ (cp r , cp^) 

is continuous, that is, outside of the inner critical layer at v = c, 
where the effects of viscosity and conductivity predominate (sec. 7a). 

To satisfy the boundary condition at the wall, cp r (, 7 i ) and 7>-|(yi) 

must be zero independently, and therefore, W = k [ T x ~ M 2 (wi - c) 2 ]= 0 
In general, c 2 ^ T x /M 2 , so that k = 0 and 


¥ = 0 , T = 0 for y _ y < 0 


(15?) 


(Cf. fig. 5.) 


By utilizing the solutions (l?3), (12^) , and (l2l+a), and equation 
(150) the di scontiriuity suffered by the Wronskian in passing from 
(y < y c ) to (y>y c ) can be calculated. In fact 1 , ' 


A¥ = 7t |cp c | 


a T _o [*^v\-j 
w* L dy V T ' Jc 


(153) 


where ® c = 9 (y c ) . Hrom (150), ¥ = kT for y = j r + 0, and therefore. 


i^c 


= it -- 


w„ 


Nl 

dy V T" j ! 


(154) 


- 1 G 


and 


a 

T = 2 k for y - y c > 0 


(155) 


In the limiting case of an incompres, 
the value calculated by Tollmien (referenc< 


idle fluid, A¥ 

10 ). 


reduces to 
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ihus^ assuming for the moment that cp G ^ 0, i f the sign of the quantity 

d /w'V" """"" 

“ \ ~~J IghJ Oglfo'-VQ} energy will pass from the me an flov to the 

disturbance; if the sign of 


d / v'\ 


»\ 1 


dy V T / 


Jc 


is negative, the mean flow 


yj-11 absorb energy from the disturban ce; 


L dy J T J i 


= 0, there 


jg-P.? ez ?.fr§P5®.. Pf energy between the mean flow and the d i sturbance . 

In the foregoing discussion, it was tacitly assumed that cp c y 0. 

By means of a proof similar to that given by Tollmien (reference 10) for 
the case of an incompressible fluid, it is not difficult to show that 
<p c cannot, vanish if cp(y) is a solution of (112) satisfying the bound- 
ary condition (JHyjJ at the wall. This presently will be done. 

, . satisfies (112), tnen cp = cpi + Ap cpp, and the be- 

havior of <p x and q> 2 in the neighborhood of the noint v = v c is 
given by (123), (124), and (124a). Since ^(j c ) ~0 and* <p 2 (y c ) 4 0 

A 2 = 0 if <p(y c ) vanishes. Now ^(y) ± s analytic and JL ( — 2— 

ay V w ~ °J 

is finite at y = y G . By direct integration of equation (112), ' 


r 


VW - (*-0) jl«* F { .. M‘J dy F liJA dy + ... I 


i 


L (w-c) 


.J T 


J 


(156) 


If yi < ,y < y 4 where y 4 is the value of ,y for which — = (w - c) 2 , 


then 


W 

T 2> 

7 ~ ”,2 - M j > 0, and the quantity in brackets in (3.56) is 


positive. Therefore, <p(y) >0 for y Q < y< y 4 and cp(y) < 0 for 

7i < y < f p(y ) can never satisfy the boundary condition <p = 0 at the 

wall, and the assumption <p(y c ) = 0 must be abandoned. 
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9. Necessary and Sufficient Conditions for the Existence of an 
Inviscid Subsonic Disturbance 


So far in the discussion of the inviscid disturbance in a compress- 
ible fluid, ,it has been assumed that solutions of equation (112) exist 
which satisfy the given types of boundary conditions. The energy criteria 
developed in section 8 not only serve to clarify the physical problem 
considerably, but also lead directly to the formulation at least of the 
necessary conditions for the existence of each of the three possible types 
of neutral inviscid disturbance . The sufficiency conditions cannot follow 
directly from energy considerations. 

This section will deal with subsonic disturbances, neutral and self- 
excited. First a necessary and sufficient condition for the existence 
of neutral subsonic disturbances is established. It is then possible 
to establish a sufficient condition for the existence of self -excited 
subsonic disturbances. However, a necessary condition has not yet been 
established. * 

(a) The Neutral Subsonic Disturbance : 

At large distances from the wall, the neutral subsonic disturb— 

8v f 1 

ance dies off like e , and 0 as y — > ® . In 

this case, no energy is transported into or out of the boundary layer 
by the disturbance, and therefore there is no net exchange of energy 
between the mean flow and the disturbance within the boundary layer 
(cf. (146 ) ) . From the results of section 8, t a C, and hence the 

quantity ~ f must vanish, if c ^ 0. If iL (JLS ) does 

not vanish for some w>l — — , the only possible neutral subsonic dis- 

M 

turbance is the one for which c = 0. When — ( — vanishes, for 

dy \ T/ 

w = Cg (say), then c equals c B for the neutral subsonic disturbance, 

d /w f \ 1 

The condition that — ( — ) must vanish for some w > 1 

dy \ T / M 

is also sufficient for the existence of a neutral subsonic disturbance. 

As in the case of an incompressible fluid, the sufficiency condition can 
be derived by means of an argument based on the fact that cp(yi; a) is 
an analytic function of a (reference l). For the purpose of this dis- 
cussion, it is convenient to deal with the disturbance equation in the 
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self adjoint form (113). Suppose no w that 


= 0 for some 


d ( w* 

. , . . dy V T . 

7 if vnere \r ~ c g (say). Then, by the necessary condition, the 
phase velocity of the neutral disturbance, if it exists^ must be equal 
to C s . Nov, . I (y) is positive continuous and bounded everywhere, and 
hence q(y) is also continuous and bounded everywhere. Equation ( 11 3 ) 
can then be integrated directly to give the relation 


7 2 2 

t (y) 9 *(y) -t (72)9* (72) ~ P ^q + -i-^cpdy 

' J 7 


(157) 


o 

By choosing the value of a 2 large enough, the quantity q + can 

always be made positive, since q(y) is bounded, low, for c = c B , 
the solution cp(y ) can be completely defined by the boundary conditions 


9 ( 72 ) 1 — c s , 9* ( 72 ) + ct \7l — M 2 (l — c p ) 2 9(y 2 ) = 0. Therefore, 

9* ( 72 ) < 0 when a > 0, and from (157), o’ (y) < 0 when a + > 0 

T 

Hence, the value of a can be chosen large enough so that 

9(yi) > 9 ( 72 ) > 0. For a = 0, however, cp(y) = y - c s and cp(y x ) < 0. 

Since 9 ( 71 ; cc.) is a bounded, continuous function of a (sec. 2), 

9(yi) must vanish for some value of <x = o, 0 > 0. For a given value of 
the Mach number, the value of c = c g is determined from the mean velocity-. 
temperature profile, and the corresponding value of the frequency a = Og 
is given by the secular equation (109). The boundary- value problem for,, 
the case of a neutral subsonic disturbance is completely solved. 


From the disturDance equation (112) and the boundary conditions, 
it can be seen that for M < 1 the singular solution cp = v, for 
c — 0, a = 0 (infinite wavelength and zero wave velocity) always 

exists, provided llm =1. If M > 1, then in 

c, a~S* o T x a 

the limiting case of infinite wavelength, (a = 0) the neutral subsonic 


disturbance becomes a neutral sonic disturbance 


is, the condition 



that 
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lim 


Wl 


' 0 M 2 J 1 — M 2 (1 - c) ; 


l 


c — - _ 
M. 


T x a 


= constant 


a o 

holds. The solution for cp is a linear 


combination of (w— c) and 


(w - c) 


f T s'! 


If M ~ 1.0, the condition 


lim „j. 
c — $> o 

a — -> o 


Wi ? c 


3 /2 


T i a 


« constant 


holds for the neutral subsonic disturbance in the limiting case of in- 
fin^te wavelength, and the singular solution <p * w exists. The sig- 
nificance of these limiting conditions will be appreciated in the inf 
vestigation of the asymptotic behavior of the a - R curve for the 
neutral subsonic disturbance in a viscous compressible fluid, which 
will be carried out in a subsequent, report. 

^ ^ lit ied and Damped Subsonic Disturbances ; 

It has been found (sec. 9) that the condition that 1 . ( Yl') 

1 dy \ t J 

vanish for some w > 1 - - is necessary and sufficient for the exist- 

ff 6 °£ a n ® u " fcral inviscid subsonic disturbance. By analogy with the 
case of an incoaro-eesme fluid, ft can be ejected that Se SSiSon 


must 


JL_ 

dy 


§■) " 0 for BOme w >:L - ~ is also sufficient for the existence 

M 


of amplified subsonic 
disturoance c = c 

s j 

value of w > 1 _ I, 

M 


disturbances ("adjacent" to the neutral subsonic 
a = a Q ). If does not vanish for some 

it appears probable that except for the disturbances 
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G ^ for M <c X. on c ~* 1 t / ?/ n x» ^ 

damped subsonic disturbance* £n er7-A 7„ ^ .' a 7 A for M - ^ onl y 

41 exist ln "t’^ 6 inviscid compressible fluid. 


io prove that the condition 


dy 


w 

T 


0 for some w 


1 

M 


18 a ; ?Et|dc3^t condition for the existence of A - 

“ r rrr, s y, 

uat da(c) = c ± i 0; (2) the sign of Ci . .“ B 

for W^ ) ™h i ^,d^f!f y ,. b8en fom4 (sec - 7) that, the boundary condition* 
terlstbo yalueTof ,X,M * '® 1?ue ~ 01 “tion between the ctanwe- 


C = c (a, M 2 ) 


(158) 


w^-,rrs f r f o on ? i ■* m:> •»«»* to & 

tw* , + ; c - 0, a - 0. In trie neighborhood of a = a 4 0 

the^e is at least one value of c for every vaii^ n-c t 1 \ * ’ 

a / c, then c 4 c ^ „ Y A mll,e of a (real), and if 

A 0 r c e> ^ jd c must ho complex; that is, 0 , 5/0 for 

the only permissible real value of c is c ' 

s 


relati A 019 n9lghi0rh °° a of » = “B>0, a and of are bnlcuely 
related, and c is an analytic function of a s Th^n n( ?\ , 

expanded in a Taylor's serie- of 1 Y - 3 ' c(a ') 1287 be 

\ o = ^2. " oi K - = a- - a 8 - around the point 


a, 


2. 

"S ' 


c = c^ + (X - \ ) ( d Y + p Y) " /d^c ^ 
^ V n-v j .T ™ ( o~ j 

4dX f ~ 4 


dX' 


s 


2 f 


(159) 


If Im f-°') 4o ny , , /d k c\ 

{ ) T V ) 01 , li Du ( )r=0V~1O9 n 1 

^ dX 4 MX C ^ ■ 9 * 3^ n ~ 1* hut 


Im 


c 


a_ c \ , 

-~rj f 0, and n is odd, then c. will always he positive for 

some value of a slightly smaller or larger than « + 1 

Of c «) » 0rt1 „x. , X J-arfaer tnan a s . ior these value: 

axis (sec! 7a ) “ “ 10n ^ y exists which is valid all along the real 
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The imaginary part of ( — ~ ^ can most readily he calculated hy 

V dX K j e 

successive differentiation of the differential equation (112) with re- 
spect to X. If <p(y) i s a characteristic function, and c and X are 


the corresponding characteristic values, 
and R" (sec.- 6) and Indeed 


d<P 

— exists in the regions R* 
dX 


dcp . ckp Sep dc 

dX ^ k 3c dX 


(160) 


For the purpose of this discussion, equation (112) can he rewritten in 
the form 


L(cp) = <p" + LL q>*_/ q + i \ ?. = 0 (l6l) 

t V T / | 

where the primes denote differentiation with respect to y. By differ 
entiating (l6l) once with respect to X, the following differential 
equation for cpy is obtained: 



When X -> X g , c — $» c g , the corresponding expressions for L(cp^) 
and L(q>) will he denoted hy I> s ( f Pxs ) and L 0 (q) g ), respectively. From 
(l60) and (l6l), 
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«S 


^8 1 ^S L 8 C<PX S ) ~ 'iPX.g^s C*P B ) J- = 


fl 


dy 


^Xs 




1 

J 




eV 


~ d 

_dc 


at- 


9 9* \ 

8 s s e 


1 AA 

.dc ^ | A 


} 


(163) 


If both sides of equation (163) are integrated between the limits y - j x 
and J - Js. along any path in the region B in the complex y— plane, an 

expression for ( — -- ^ is obtained. Consider first the integral of the 
v dX 'b 

left-hand side of equation (163): 



9 f <Pa 
s Xs 


\ 


Js 


' 1 yi 


( 1 6 k) 


Since 9 ( y j.; c, X) = 0 is an identity in X, 
and the integrated expression vanishes at y^. 


9 8 (.7x) = 9x s (.7i) = 0, 
At the upper limit. 


9*s C, X) - - a s A -If (1 - c s r q>g(y a ; c B , X B ), 


9 f x s (jz, c B , X g ) = - a s -v/l - M (l - c s ) s cp ?kS (y 2 ; c s , X g ) 


” 9g 


(y s ; c«i, X e ) f a Jl - m" (l - cf 1 


v. dX 


JB 
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Substituting these relations into (l6b) finally results in 


1 e ( Wb - %) 


I , / do \ 

— 8*1 + 83 f — - j 

* ^ dX 's 

yi 


(165) 


where, ax and a, 2 are real constants. 

The integral of the right-hand side of equation (l6M is 


,^2 ffl 


T 


dy + 


7x 


dc^j 

d */s 


I 72 


9, 


S/M 


m + ( 

dc y 3 ' 


7i 


q + .M) | *3 (Ml 

S T ' ** s ^dc i ^ 


9 S 9 £ 


f d V. 




dc l 


s - 


tPq"* /dC\ 

ay = / JL_ dy + f ) I x 


y 1 


VdXA 


(166) 


Equations (l64), (l65),.and (l 66) yield 


dcN 

dxv c 


Jf2 

2 

9s dy + a x 

'7i t 

Ii — a? 


(167) 


72 9s 2 


In evaluating the integrals Ii and / dy, integrate along the real 


axis, except for the term 


>72 


dy. Indeed, all the other 

v dc y s 


'7i 
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integrals have real and finite integrands along the real axis. Thus, 
the imaginary term in (l6 r f) can occur only with the Integral 


pys 2 

/ — dy, the integrand of which becomes infinite at y = y . 

J Sic y B ° 

J y i 

By expanding the integrand in power series in the neighborhood of y = y c , 
there is obtained 


Im ( I x 

L 




= ImJ 


sc +a Pec < p , s C (y-y c J + •**)_( ^y c ) + * 

Jyi I.’ef (y-y c ) a f 1+ (y-y,,)+ . . . V 

C ■ ° L 2w» - j 


1 

-1 dy 


J 

(168) 


from which. 


Im (l 2 ^ 


<p, 


c 


( ^s w ‘ ) "c 

( v f .*) 2 


(169) 


If 9_ (y ) f s a characteristic function, cp can never vanish (sec. 8), 

s BC 

ImY^-^ 0, provided ■(! s w t )" ^ 0; and 

c i must be positive for some value of a slightly smal ler o r larger 

than a a , 

The restriction that (tcW')" must not vanish can easily be re- 

0 c 

moved by an extension of the foregoing argument. By the physical nature 
of the mean velocity- temperature profile, if (I w f )* r = 0 the quantity 

1 gW’ must have a true extremum at the point y = y c w » c s and not an 


Hence, Im 


r) 


i 0, 
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inflection; in other words, (IgW 1 )* must have a zero of odd order at 
the point y = y c . Therefore, 



nm+x 

( y-y c ) 

(2m+l) * 



♦ } SI » 0 y 1 } 2, • . . 


( 170 ) 


do 

From ( 169 ), it can he seen that Im ( — ^ = 0 for m ^ } . The differ— 

V ye ' "■ 

ential equation (162) for cp^ is regular in the vicinity of the point 
y = y c , and therefore cp> v (y) is real. By differentiating (162) suc- 
cessively with respect to X, differential equations are obtained for 

which are all regular in the vicinity of the point 
Consequently, cpy ^ is real for k < 2m+l. In 


d^q) £3^ 

dX 2 ' d)J : ’ ’ ' dXk' 

y = y C if k < 2m+l. 


k 


the expression for ( !L£ 


dX k / s 


a term of the form 


k! (^f 

W / 


r 12 t s 2 


(7T77 c+1 ^ b8 " ) 


(ie* 




dj 


ji 


always appears by analogy with (168). (All the other terms in the ex- 
/d k c\ 

pression ior f — p) are always real.) By virtue of (170), the im- 

&V y s 

aginary part of this term vanishes if k < 2 m + 1. However, if k = Pm-1 
the imaginary part of this term is " ' 


n 


, am+x f I 
dc x !<p c ! 


dX 


(w» ) 2BH2 v 

c 


/ \ ( am+ 2 ) 

(is*') i 0 

y c 


2m+x 

Therefore, Im ( -y ^ 0, and Cj_> 0 


dX 


am+i 


for some value of a slightly 


larger or smaller than a g . The proof of the existence of amplified 
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subsonic disturbances adjacent to the neutral subsonic disturbance 
c = c g , a = a g is thus complete. 


It is quite difficult to give a rigorous proof of the existence of 
amplified, subsonic disturbances adjacent to the neutral disturbance c = 0 
a ~ chiefly because c is not an analytic function of a in the 
vicinity of this point. (See secs. 7 and 7c.) Although it does not seem 
worth while to discuss the details, it can be shown by a method similar 


to that utilized by Tollmien (reference 10) that ( > 

^ OX ' 


is real, 


c ~ o 


Im 


d 2 c 

dX 2 


is positive if 


c=o 


' a 

d y 


T J ' 


> 0, and 


Ji 


dX 


2 'c=o 


is unfortunately logarithmically infinite. The argument in this case 
is therefore inconclusive. Of course, from the asymptotic behavior of 
the neutral a-R curve, it should be possible to see that amplified 
subsonic disturbances do actually exist in the neighborhood of the neu- 
ral disturbance c = 0, a = 0, if the neutral disturbance c = c<?, 
a = a 8 exists. 


10. Some Further Discussions of Inviscid Disturbances 

So far, only subsonic disturbances which are neutral or nearly neu- 

ne iffhWh nS US 86 ?! 4 ^ hese disturbances correspond to the immediate 
g . g l g hb o rhood of the positive real axis of the complex Q -plane (fig.“ 3 'J 7 " 

iu has not yet been possible to get any result regarding general non- 

modes, except that they possess the property of being either self. 
a J d out f°ing, or damped and incoming, and that no sharp change 
in property would be expected in passing from subsonic to sonic and to 
S disturbances. The neutral sonic and supersonic disturbances, 
however, do enjoy a special position. The former corresponds to the 

branch point of J q at the origin, and the latter corresponds to the 
cut drawn in the Q -plane to separate the two solutions exp (± a */?!>). 


In the following sections, the neutral 


. considered. The transfer of energy is made the basis of this 

ves gation. In section 10c, the case of the neutral sonic disturb- 
ance is discussed briefly. 


Of course any other cut might have been used. This particular 
one, however, has the desirable property that one of the solutions 

exp i: a y 1> is in general ruled out by physical requirements. 
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(a. ) Nec essary Conditions f o r th e E xistence of a Heutral Supe r s 
Disturbance : 


onic 


The results of the investigation of the energy "balance for any type 
o neutral inviscid disturbance lead directly to the necessary conditions 
ior the existence of the neutral supersonic disturbance: 

(l) If the mean flow in the boundary layer absorbs energy from the 


disturbance. 


A (?!') 

dy V T J 


j y=y c 


0, c > c„. The amplitude of the re- 


flected wave must be less than the amplitude of the incident wave. 

, ^ ^ ere no exchange of energy between the mean flow and 

the disturbance, c = c s . The amplitudes of the incident and reflected 
waves are equal. 

(3) If the disturbance absorbs energy from the mean flow in the 

r 


l dy x T 


boundary layer, j A ( > °> c < c q* Ths amplitude of the 

\ m / : 


y=y c 


reflected wave must be 

Of course, c < 1 _A 

M* 


greater than the amplitude of the incident wave, 
and a is arbitrary. The stationary Mach waves 


c = 0 (a arbitrary) can always exist. 


The necessary conditions for the existence of a pure outgoing or a 
pure incoming wave will be discussed in connection with the reflection 
and absorption of the neutral supersonic disturbance. (See sec. 10b.) 
-Formulation of the euj^icj^nt conditions in this special case has proved 
to be a formidable task. In general, q(y) is not bounded at the "point 
w - c, and l (y) is not bounded at the point T - M 2 (w ~ c) s . Con- 
sequently, it is difficult to determine the sign of c?(y x ; a) for large 
values of a, ^and it has not yet been possible to carry through the type 

of argument which served in the case of a subsonic disturbance. (See 
sec. 9a.) 


^ Deflection an d Absorption of the Neutral Supersonic Disturb ance: 

By the action of the viscous forces within the critical layer at 
w - c (c f c g )} e relative phase shift is produced between u* f and 
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* 

V *S and the shear stress t = — — — increases (or decreases) rap— 

Po*(u 0 *) a 

idly from zero for (y _ y ) < 0 to the value - k for (y - y„ ) > 0 

2 1 c 

(secs. 4 and 6). Thus, the critical layer takes the place, in a sense, 
of a wavy wall or irregular solid boundary, in reinforcing or partially 

canceling ("depending on the sign of I .5 ( "'ll ^ 1 ^ an incoming 
" !_ dy 'T / J c / 

disturbance during the process of reflection. 

From equation (l48), the shear stress t = ^ k is also equal to 

, , . a 1 . icy . — icoy 

the expression — _3m(cp»<pJ. Since <p - A© + Be 

2 T - M s (w - c)^ 

for y >> 1 (sec. T) 


T = °L k 
2 


2 J M 2 (j 


c) a - 1 


B j _ A 


= 1 


(171) 


and hence, from equation (154), 


J M 2 (l -of 


fbfkn .* ,<p i 




f L 


dy T - 


(172) 


By making use of equation (172) and the additional relation 


19 f ( 72 ) 


+ C0 s j Cp ( V 2 ) 


r 


= 2or ' ' | A 1 + B 


n 


Bi 


(173) 


an expression is obtained for the "reflectivity” K - defined as 


the ratio of the energy carried out of the boundary layer by the reflected 
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wave to the energy brought into the boundary layer by the incident wave. 
It is found that 


where 


i - j 


(m) 


It follows that 


2jto- 


_ a 


1 d 
•v 1 dy 


1 

V a 


7 = Jc 


9'(ya) 


9c 


4- ur 


9(ya) 


% 


K = 1 when ,j 


> 


0 


( I J T< 1) 


(175) 


The necessary condition for the existence of a pure outgoing or incoming 
wave is 



(176) 


In this case, the boundary condition at y - y 2 is 


and 


cp'(ys) ± da*p(y 2 ) = 0 


«oj 

a2 


j 1 d f w f 

L w* dy ^ T 



1 9 (. 72 ) 


2 




( 177 ) 
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In section 7c, it was remarked that a solution of the disturbance 
equation (112) satisfying the boundary condition cp = 0 at the wall 
(7 = Jit or y = 0) could always be found for arbitrary values of c 
and a in the case of the supersonic disturbance. Such a solution is 
(cf. (131) 


so that 


00 

<P<y) ■* (w - c) y a 2n k 2n +i (yj c, M 2 ) 
n=o 


(178) 


Vfos) = (l - c) y a 2n k 2n +i (jsi c, M 2 ) ( 179 ) 

n=o 


^'(ys) = °) y a 2n k* 2 n+x (ya; c, M 2 ) 

n-o 


( 180 ) 


9c = <p(yc) 



lim 

y^y c 


(w-c) 2 y a 2n k’sn+i (y c i c, M 2 ) ] 
n~o 





(w-c) 2 

T 


dy 




r T 

\(w~c)‘ 


- M 2 1 

J 


dy 


1 

•i (l8l) 
/ 


It is assumed that the mean velocity temperature profile in any 
particular case is known for each Mach number. Subject only to the re- 
striction c < 1 - i, the reflectivity can be calculated for a series of 

M - 

suitable real values of c. These calculations should give some 
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indication of the conditions under which a pure outgoing wave or a pure 
incoming wave can exist. 


Necessary and Sufficient Conditions for the Existence of a 
Neutral Sonic Disturbance : 

If the physical condition that both <p and *r/p must be bounded 
as y — 5> co ie imposed, then in the case of the neutral sonic disturb- 
ance, <p and <p* — > 0 very rapidly as y «. (sec. 7c) and no 
energy can be transported into or out of the boundary layer by the dis- 

turbance. The necessary condition for the existence of a neutral sonic 
disturbance Is therefore (sec. 8) 



(a ^ 0 ) 


(132) 


j j case ^Ne neutral subsonic disturbance (sec. 9a), the 
condition c = c s is not entirely sufficient for the existence of a 
neutral sonic disturbance. Because the physical significance of this 
sonic disturbance is not yet clear, it does not seem worth while to 
discuss this problem in great detail, although some mathematical results 
have been obtained. A brief sketch of the arguments and results will be 
given here. From equations (117) to (120), and equation (138), the so- 
lution of the differential equation (112) for a = 0 which has the cor- 
rect asymptotic behavior in this case must be 


Therefore, 


9(y) = (w _ 


c) 



T 

(7T7T 2 " 


- M‘ 


a 

j 


dy 


(183) 


<p(0) 



a = 0 


(184) 
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and hence, 


according as 


<p(o) 5 0 



f l 

\ (w - c) 2 


- M" 


dy < 0 


(165) 


011 the otiier hand > 8X1 argument almost identical with that utilized 
in the case of the subsonic disturbance (sec. 9a) shows that <p(0; a) > 0 
or large values of a, if c - c 0 . Since cp(0^ a) is n bounded con- 
tinuous function of a, if 



( - T _ M 2 1 

l (w - c) 2 J 


dy > 0 


cp(0, a) must vanish for some value of a, > 0; 


if 


r°° f 1 M 2 1 dy = 0 

J l(w - c) 2 S 
^0 


a non-trivial solution exists for a = 0; if 


r * 00 




_ M' 


.2 1 


dy < 


0 


it must first be determined whether or not cp(0, a) > cp(0 0) 
a, before any definite conclusion can be drawn. By employing 
tion of the oscillation, or comparison theorem (reference ll) 
shown that <p(0, a 2 ) > <p(0, a x ), if a 2 >a x , and therefore 
is a monotonic increasing function of a. Hence, if 


for all 
a modifies- 


it can be 
9(0, a) 
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00 



f 1 

> (w- c) a 


- M 2 


| ay< o 


no solution of this type exists 1 , with a 4 0 . (See end of sec. 9a.) 


11. Concluding Discussions 

The above investigations of energy relations and the necessary and 
sufficient conditions for the existence of certain types of disturb- 
ance, though incomplete, serve to give a general understanding of the 
stability problem in an inviseid fluid. Before proceeding to include the 
effect of viscosity, the significance of the results will be discussed 
somewhat in detail, 

The distribution of the density of angular momentum across the 

_* du* 

boundary layer is unstable if the quantity p has an extremum 

dy * 

for some positive value of u* > u* ^1 — i ^ , where M is the Mach 

number for the mean flow outside the boundary layer. From the equations 
of mean motion, it is not difficult to show that the quantity 


1 (pv<) = 

dy dy V T 


will always vanish for some value of w > 0 if the solid boundary is 
insulated, or if heat is being transferred to the fluid across the solid 
boundary. This can be seen as follows: 


The dynamical equation for the mean motion 


_* / — * CiU — * du 

p ( U + V 

' c)X* 3y* 


dy* 


_* 



X 


Subsequent investigation has shown that this is the case. 
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gives 


O U N 

m 


. du 

dy ' 2 ' 

) rr& 

1 T 

My* A 

dy* 


if 


M 
— * 


— *X- 

T 

J 

To 


(m = O .76 for air) 


where the subscript " 1 " refers to the wall. If heat is transferred to 

—'-'a - — y. 

the fluid at the wall, ( zLS\ < 0 and hence ( SDLA > 0 . There- 
fore, the quantity 


— * 

1 du 


¥r _ — -X- *• 

1 / B 2 u du 1 BT 


By* \ T* By* / V By* 2 By* T* By* 


must he positive for y* = yi*. Thus, 


1 du 


T dy* 


increases from some posi- 


tive value as y increases from y x *. But also it is known that it 

approaches zero as y* becomes infinite. Hence, ~ T BL has a maximum 

T ' by* 

at some point y*> yx*; that is, -!Lf vanishes for some w > 0. 

dy V T ) 


If the solid boundary is insulated. 


dl 1 " 

dy* 


/d 2 u \ 
= 0, and ( ) 

a ^ oy* E ' 


“0. The above argument yields no conclusive result. By differentiating 
the dynamical equation once more with respect to y*, and utilizing the 
equation of continuity for mean motion, 
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dx* by* 


U* + y“ Sp 


ax* 


by 


it is found that 


= om /du*V n 

\ X- 3/ .T'X- V v j 

'by ' x CpTi \3y V 


d 2 "* 

Thus, must - De positive for some y* > y\y and since < o, 

by* 

for y* < y*, it follows that the quantity ^ Q. |!) muS t be posi- 

case L°L 8 ^ ?"* ' Hence * the essential conditions for the last 

se also hold m this case. The same conclusion is therefore obtained. 

However, if heat ^withdrawn from the fluid at the solid boundary 
/ 


dT x . 

j > 0 and 

by * / 1 


d 2 u* \ ^ 

< 0. The signs of the quantities 


and 


> 2 ~* 
a u 


by 


by* 


^*2 wili remain unchanged as y* increases from y* x in general. 3 
Hence, the quantity remains negative and will not vanish. 

is lu T J? < 1# if the boundar 7 is insulated or if heat 

in ^ nof ht int ° th ® fluid > dt is certain that the laminar boundary layer 

Blasiu® boundiv 1 °i flUid - Wi11 ^ relativel y less stable than the isothermal 
d ^ la ?' 6r in an ^compressible fluid, as far as the inertial 

more s tab L^Al though ll iB , teken from the fluid > the flow will be 
more stable. Although these conclusions can probably be extended to the 


Except when the quantity 


(< - 5) o 1 ' 


TO ~ Tl ) ]/ ^ l0 V6ry ™ Sl1 ’ Where T ! 
perature . 


-vf 

i 

s 


stagnation tem- 
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case when the velocity of the mean flow outside the boundary layer is 
only slightly supersonic, no statement can as yet be made for the general 
supersonic case. 1 

However, the critical Reynolds number defined in terms of free— stream 
quantities may not necessarily be decreased by heating the solid boundary. 
For in the viscous solutions of section 3, it is the kinematic coeffi- 
cient of viscosity near the solid boundary that enters. This coefficient 
is increased by heating, thus leading to greater stability. Whether the 
minimum critical Reynolds number for any compressible— fluid boundary layer 
at any Mach number will be greater or less than the value for the Blasius 
profile can be determined only by actual calculation. This question will 
be settled for several representative cases in a forthcoming report by 
some numerical work following methods to be discussed in the next part 
of this report. 

a^recent report (reference 6), Allen and Nitzberg suggested 
that the proper" Reynolds number should be based upon the kinematic 
viscosity at the solid boundary. For small values of c*, this 
is not very different from that at the critical layer. However, 

— £- „ * 

they have assumed that the critical Reynolds number \ is 

V v 1 1 /cr 

equal to the critical Reynolds number for the Blasius pi-ofile. For the 
case of insulated solid boundaries (e. g,, airfoil surfaces), their 

value of i^u 6 /{? 3 Q may therefore be too high. 


i — ~ ‘ 

For example, for M > 1.5* it may not be possible for a subsonic 
characteristic-oscillation to exist in certain cases, because in addition 
to satisfying the equations of motion and the boundary conditions, it 
must also satisfy the condition c*> 1- l/M. 
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III - STABILITY IN A "VISCOUS CONDUCTIVE GAS 
X'd. General Considerations and Methods of Numerical 
Calculations for the Stability in a Viscous Fluid 


The foregoing inviscid investigations serve to illustrate the gener- 
al behavior of the pressure and inertial forces in the control of the 
stability of the flow of a compressible fluid. These results can there- 
fore be used as a guide in the investigation of the stability in a real 
fluid at large Eeynolds numbers. In the case of , the incompressible 
fluid, very valuable information has been obtained by consideration of 
a iA9Al£ ig , at i on of the results in an inviscid fluid by the effect of vis- 
cosity.. The general conclusion has been reached that the effect of vis- 
cosity is essentially destabilizing at very large Eeynolds numbers; and 
it has been possible to obtain the asymptotic behavior of the neutral 
stability curve for large values of the Eeynolds number, also to give a 
Quick approximate estimation of the minimum critical Eeynolds number and 
indeed to compute the complete curve of neutral stability. In the present 
case, corresponding developments should also be possible", but the results 
evidently depend upon the Mach number. Any computation of the curve of 
neutral stability must be carried out for each value of the Mach number of 
the free stream. 

Owing to the limitations of time, it has not been possible to carry 
out these computations. The authors, however, laid down the general plan 
of the calculation of the neutral curve of stability, and repeated the 
calculation of Tietjens function. Some of the numerical values turn out 
to be slightly different from those originally given by Tietjene (table I, 
fig. 6). They agree very closely with the results of Schlichting's later 
calculations (table 2 9 p. 73 > roforonco 5 ) ♦ 

A method of numerical calculation very similar to that used in the 
incompressible case will be outlined below. It enables the curve of 
neutral stability to be computed for each Mach number as soon as the dis- 
tributions of velocity and temperature are known for that Mach number. 
Several such distributions have been obtained by Karman and Tsien 
(reference 12), by Crocco (reference 13), by Emmons and Brainerd (refer- 
ence lk), and by Hantzsche and Wendt (reference 15). 

Mgbhgd, pf numerical calculation .- The calculation of the neutral 
curve depends upon a proper evaluation of the function E(a, c,M^} 
occurring in (10 5). According to (108) , its evaluation den ends upon 
the evaluation of cp ± j and cp^ (i, j = 1,2) . To evaluate these 

functions, the inviscid solutions (117) and (118) are used. After a 
little calculation, there is obtained 
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Substituting ( 187 ) into (l08) gives 

w[ (Cp22 + 0 ?22) 
E(a,c,M 2 ) = 


f!ii r 

wl {<pJ, 2 + ptp a2 ) + 1 i l <fl s + P¥ 1; , 


CJ 

where X = X(c ) is defined by 

vi (ji - y c ) = -c (l + X) 
By introducing the function 

F(z) = 


+ X 


(189) 


1 - F(z) 

and using the relation ( 189 )# ( 105 ) may be reduced to the form 


( 190 ) 


(191) 


F(z) = 


( l + X) (n + il v) 
1 + X (u + iv) 


with u + iy = 1 + 


« 1 \ 

V x C (cp22 *f ]3Cp2 o) 


( 192 ) 


It is noted that the quantities T 22 and q> 22 involve the Integrals 
K 2 n and K 2n+i , These integrals involve 

ya 


/ 


r t 


■-^yi l -(w - o)* 1 i 


M®1 dy = K x 




f T 

•(w. - c ) ' 


*1 


M~> dy 


in the first step of integration. By substituting the, right-hand-side 
expression into cp 22 + pcp 22 , it is not difficult to verify that the 
terras involving Kp combine to give cpi 2 + 09i2* Hence, it is conven- 
ient to write ■ - 

922 + 0922 = (9X2 + 39l2) Hi -’V, 

Substituting this into (192) gives finally (with y 2 - yi = l) 


^ 2 


wi c r 

u + iv = 1 + ■— — < 


T] 


T ,1 

Fvpp ^ - M j 


Wi c i> 


(Ti)(9i2 + P9x2) 


(193) 


y.2 


yx 


n 


The function 


• dy and the integrals involved in 9xs> 9 x 2 # 


■J (w - c) 2 , 

lx 

and $ may be evaluated by methods similar to tdics© used in the 
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incompressible case. The significance of bringing the final equation in 
the form (193) is that the imaginary part of the right-hand side 


7s 

is mainly contributed by the term involving J* T(w - c)" 2 dy the 

' 7i 

imaginary part of which is - S . 5si. This can be easily calcu- 

* ~ w ! T* f 


w 1 T 
c x c 


lated. Thus, using the fact that X is usually very small gives, 
approximately, 

Fi(z) « -»t Is - S\ V i_£ 

*'o\ *'c T J T i 


( 19 ^) 


where Fj.(z) is the imaginary part of F(z). The relation (194) would 
give a correspondence between c and z. From this, the value of <xR 
can be eaeily calculated by means of (107). For more accurate calcula- 
tions, use is made of the relations 


Fr(z) 

Ji(z) 


(l + X) ^u(l + Xu) - Xv 2 j |(l + Xu) 2 + (Xv) 2 | 


(l + X)v 4(1 + Xu)"" + (Xv) 2 


/ 


(195) 


J 


where F r (?) is the real part of J(z) . Using (194) as the initial ap 

proximation, a method of successive approximations can be developed 
exactly as in the incompressible case for the calculation of a, R 
for given values of e. The complete numerical calculation will be 
carried out for a few typical cases. 


California Institute of Technology, 

Pasadena, Calif., November 15, 1944. 
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Table I.- Functions F(z) and F(z) 


— 

z 

F 

x r 

F i 

F 

±r 

l i 

1.0 

0.89161 

-0.35025 

0.80630 

•-2.60557 

1.2 

.78969 

-.27310 

1.77012 

-2.29854 

1.4 

.71970 

-.21213 

2.26836 

-1.71669 

1.6 

.66931 

-.16009 

2.44985 

-1.18600 

1.8 

.63143 

-.11274 

2.48104 

-.75892 

2.0 

.60144 

-.06741 

• 2.43927 

-.41253 

2.2 

•57599 

-.02226 

2.35196 

-.12348 

2.4 

.55230 

-.02395 

2.22724 

+.11916 

2.6 

.52773 

-.07203 

2.06929 

.31558 

2.8 

. 49952 

. 

+.12220 

1.88566 

.46043 

3.0 

. 46456 

.17391 

1.68938 

.54872 

3.2 

.41947 

.22520 

1.49726 

.58082 

3.4 

.36110 

.27193 

1.32516 

.56401 

3.6 

.28802 

.30705 

1.18429 

.51074 

3.8 

.20352 

.32130 

1.07982 

.43560 

4.0 

.11800 

.30721 

1.01118 

.35220 

4.2 

.04698 

.26559 

.97361 

.27133 

4.4 

.00240 

.20811 

.96056 

.20038 

4.6 

.02160 

.14475 

.95989 

.13601 

4.8 

.01477 

.09875 

.97659 

.09503 
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